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BEAM EFFLUX FROM A 30 CENTIMETER THRUSTER 
1 . 0 INTRODUCTION 

The operation of a mercury electron bombardment ion thruster 

results in the generation and acceleration of a beam of energetic thrust 

ions and accompanying neutralizing electrons moving along the thrusting 

axis in a generally well directed flow. The flux patterns of these thrust 
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ions have been determined in a large number of experimental programs 
(^Q^ 2 .ng the development of the mercury bombardment thruster. In addition 
to the thrust ion currents, a series of particle effluxes from the 
thruster emerge with varying particle release rates, cones of divergence, 
are energies. The released species include both neutral and ionized 
mercury, and, in smaller amounts, both neutral and ionized metal atoms. 

,8 

To distinguish the various mercury ion species, Staggs, et al 
introduced the notation of Group I, Group II, Group III, and Group IV ions. 
Group I ions, described above, are energetic ions which have been accelerated 
by the total potential difference between the bombardment discharge and 
the neutralized thrust beam and have trajectories which are generally 
contained within a cone of divergence from the thrust axis with 30° 
half-angle. Smaller quantities of Group I ions occur at larger divergence 
angles. The large angular divergence regime also contains the ion species 
designated Group II and IV. Group II ions are the result of charge transfer 
processes between mercury atoms and ions in the region between the 
bombardment discharge and the accelerator grid and, further, in those 
portions of this interspace in which the potential is positive with 
respect to the neutralized thrust beam plasma potential, ^he ions formed 
in this charge transfer process, thus, are capable of acceleration and 
release into the thrust beam, although, as will be seen, their trajectories 
are more broadly distributed in space as a result of non-optimized ion 
optics and generally high decel-accel ratios. Group III ions are the 
result of ion-atom charge transfer reactions in the bombardment discharge- 
to-tlirust beam interspace regions in which potential is negative with 
respec . to the thrust beam potential. Since escape into the thrust beam 
is energetically forbidden, the ions are collected at the accelerator 
and are not observable to thrust beam diagnostic probes. The remaining 
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ion species, Group IV, results from charge transfer reactions between 
thrust ions and mercury atoms escaping from the bombardment discharge. 

The -initial energy of these ions is, essentially, the thermal energy 
of the escaping neutral atom, and the resultant ion is subsequently 
acted on by the weak internal electric field structure in the neutralized 
plasma thrust beam. Because these electric fields are widely divergent, 
the Group IV ions emerge over a broad cone of directions. 

The interest of this beam efflux measurements program is the 

quantity and angular divergence patterns of Groups I, II, and IV ions, 

for the divergence angle regime from 0® to greater than 90®. The program 

here continues and extends earlier measurements of these ion flux patterns 
1 3 4 

with both 20 cm and 30 cm * diameter mercury engines. The engine body 
utilized in this program is the same as that utilized earlier in the 30 cm 
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beam measurements . The accelerator grids, however, have been changed 
to produce a more collimated thrust ion flow. 

This report will describe the experimental facilities including 
the testing chamber, collectors, shrouds, and engine diagnostic array, 
and the thruster in Section 2.0. Section 3.0 will review "facility 
effects" which are present as a result of thruster operation in a bounded 
(laboratory) geometry and will include analyses of the form and possible 
extent of these facilities effects. Section 4.0 will present low energy 
(Group IV) ion flux measurements, while Section 5.0 will describe models 
and analyses of this low energy ion plume. Energetic ion measurements 
are described in Section 6.0, with a discussion of models and analyses 
of these ions given in Section 7.0. The possible use of low energy ion 
flux measurements as an in-flight diagnosis of thrust performance is 
discussed in Section 8.0. This assessment of Group IV ion flux density 
as a thruster diagnostic will also consider the effects of Group II ions 
whose presence at high angles creates a "noise" signal to the Group IV 
determination, and whose combined presence and energy creates potential 
problems in diagnostic probe erosion and secondary material transport 
and deposition. 


2 


2.0 EXPERIMENTAL FACILITIES 

The testing chamber and diagnostic probes used in the thruster 
plume measurements are illustrated in Figures 1 through 8. Figure 1 
shows thruster placement in the 5’ x 11' chamber, and the location and 
size of the upper and lower shrouds and the beam collector. The shrouds 
and collector are electrically isolated from each other and from chamber 
ground. This separate electrical isolation results in a thrust beam 
neutralization condition in which overall electron flow from the neutralizer 
must equal overall ion current from the thruster and in which electron 
and ion currents must balance at each shroud and at the collector. The 
electric field pattern in the neutralized thrust beam plasma, thus, is 
an accurate simulation of engine operation conditions in space. 

Chamber pressure in the 5’ x 11’ facility depends upon the 
level of thrust ion current and upon the level of refrigerations in the 
shrouds and collector. For maximum liquid nitrogen cooling of the^shrouds 
and collectors, chamber pressure remains in the range from 2 x 10 Torr 
to 5 X 10"^ Torr. The effects of this ambient chamber pressure will be 
discussed more thoroughly in the sections dealing with facility effects 
and in the low energy ion plume measurements and modeling. 

Figure 2 provides a view, in isometric, of the thruster, 
and the diagnostic probe array. As an additional aid in 
the visualization of the probe array, Figure 3 illustrates the location 
of the probe mounting shafts, viewed along the axis of the thruster 
and the test chamber. The probe mounting shaft location and the specific 
method of coupling of each probe to its mounting shaft will determine 
the available range of probe location relative to the thruster and 
principal directions in which the probe scans the thruster plume. 

Figure 4 provides details of the Engine and 1-1/2" probes. 

The Engine is a two element Faraday cup. The outer case and grid are 
electrically connected and are biased, generally, at a potential negative 
with respect to the thrust beam plasma to prevent electrons from that 
plasma from moving to the collector and to suppress secondary electrons 
emitted from the collector by impact of the thrust ions. The rotation 
of the Engine mounting shaft causes the Engine probe to swing on 
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Figure 3. View Along Thrust Beam (z) Axis Illustrating Radial Distance 
and Azimuthal Position of Diagnostic Probe Mounting Rods. Probe Motion 
and Position Indicated for Engine J , 1-1/2" J , and Swinging J . 

Probe Radial Position for 4" J , Piggyback J. , and J Weasel I and 
Weasel II also Indicated. ^ 
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Configuration on J Weasel I and II Probes. 



an arc which passes through the axis of the thrust beam, thus allowing 
the probe to determine thrust ion current density as a function of 
r and z in a cylindrical coordinate system, (r, z, '(>). whose z axis is the 
beam axis. When the case and the collector are connected to each other 
and isolated from chamber ground by a high impedance (10 megohms), the 
probe acts as a floating probe to determine the floating potential in the 
neutralized thrust beam for the same r and z locations for which the 
probe, used as a Faraday cup, determined ion current density. 

The 1-1/2" probe is also illustrated in Figure 4. This 
three element probe consists of an outer grid and case, an inner grid, 
and a collector and, by setting a negative bias on the outer grid and 
varying the bias on the inner grid, permits the use of the probe as 
a retarding potential analyzer. The probe, thus, not only determines 
ion current density, but, through retarding potential analysis permits 
a determination of the energy spectrum of arriving ions. The mounting 
and motion of this probe is the same as that used in the Engine 
except that a longer mounting shaft is used, thus permitting ion current 
density measurements to be carried out over a larger interval in axial 
distance, z, than is possible with Engine J_^. 

Figure 5 illustrates the motion and construction details of 
the Swinging-J^ probe. The axis of this four element probe (outer case 
and aperture, inner and outer grids, and collector) intersects the face 
of the thruster at the thruster axis (r = z = 0) and remains fixed on that 
point as the probe is rotated. The probe motion, thus, is in 6, where 
0 is the angle of divergence from the thrust beam axis, and results 
from a drive shaft through the Swinging probe mount and appropriate 
gears and coupling to cause the probe arm and probe motion. The function 
of the probe is ion current density as a function of divergence angle 
and ion energy. 

Figures 6 and 7 provide views of the 4" probe and 
Piggyback probe package. Figure 6 (see also Figure 2) shows the 
probe package as viewed from various directions. Details of the collector 
inner grid, and outer grid (and case) of the three element 4" and the 
collector, inner grid, middle grid, outer grid (and aperture and case) 
of the four element Piggyback are given in Figure 7. The probes, 
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used as retarding potential analyzers, provide a measure of ion current 
density as functions of probe position and ion energy. The 4" has 
a comparatively broad acceptance cone of directions and in its usual 
orientation determines the totality of ion current impinging on a surface 
whose surface normal is in the radial direction (in the cylindrical 
coordinate system described earlier) and which intersects the thrust 
beam axis at location z. Through motion of the 4" mounting shaft 
the location in z is varied. The face of the 4" is at r 28 cm. 

Because of the width of the collector plate relative to the 
spacing between grids and because of collector placement within the 
outer case, the 4" determines the total current of ions moving in 
the 0 range from 'v 0® to 180® and intercepting the cylindrical surface, 

r ® ’^probe* ^ “ ^probe’ angular range of the Piggyback 

probe is more restricted, however, and it serves primarily to measure 
ion currents in the "backward" (6 180®) direction, again, as a function 

of probe location z. The principal interest in these backward streaming 
ions is in terms of facilities effects, since operation of the thruster 
in space does not, in general, lead to any significant levels of ion 
current in this reverse direction. 

A final figure illustrating probes in the array is Figure 8, 
which shows details of the Weasel 1 and II probes. The Weasel 
contains 5 separate collectors. The collectors lie behind an inner and 
outer grid, each of which is separately biasable, and inside an outer 
case and grid which is also separately biasable. The probe moves in the 
2 direction through the motion of the mounting shaft. In common usage 
the axes of the various collectors intersect the thruster axis. In this 
orientation the probe determines ion current density as a function of 
ion direction of travel in 0. The angular width of each collector (in 
conjunction with the aperture in the case) is ■ 30®. By exercising 
probe motion and bias it is possible, In principle, to determine ion 
current density as a function of 0 and ion energy at position z for z 
along a cylinder of radius r , . 
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3.0 FACILITIES EFFECTS 

3 . 1 gen eral Consldorations 

TUe operation of an ion thruster results in the release of 
several forms of particles, including both charged and neutral species 
with energies varying over a broad range for both forms of particles, 
and a variety of RF and optical emissions. The measurement of these 
wave and particle fluxes constitutes the principal tasks in determining 
interactive effects between the thruster, the ambient space, and a 
spacecraft. Fcr laboratory determinations of interactive effects, 
measurements (of necessity) are conducted in the bounded geometry of 
the laboratory test chamber, and the term "facilities effects" can be 
interpreted to mean any alteration of a measurement of thruster wave 
and particle release resulting from the presence of the testing chamber. 

In this program report, the concern of facilities effects will be more 
narrowly limited to measurements of the charged species. 

In principle it should be possible, by increases in chamber 
size or by reductions in chamber background pressure, to lower the level 
of facility effect particle fluxes below the levels of "genuine" particle 
emissions from the thruster. In practice, and depending upon the specific 
details of a given measurement, the reduction of facilities effects to 
negligible levels may not be easy, and, in some instances, may not be 
possible in any real and practical sense. Key elements in assessing 
possible influences of facility effects is the spatial location of a 
field point in question, the direction of arrival of particles examined 
by a probe, and the energies of the particles examined. Measurements of 
thrust ions for probes situated on or near the thruster axis should not 
be influenced significantly by facility presence. Measurements of very 
dilute fluxes of charged particles at either large spatial separation 
from the primary beam or at very high angles of divergence from the 
thrust beam axis may be significantly influenced or even dominated by 
facilities effects. The discussion in this and the following sections 
will attempt to identify some of the reactions and locations for which 
facility effects are of concern. 
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3.2 Neutral Particle Species 

In the period before onset of thruster operation, neutral 
particle species include the common gases for laboratory chambers (N^, 

O 2 , H 2 O, and hydrocarbons (these last from diffusion pump backstreaming , 
vacuum grease coverings of seals) and umpumped mercury atoms). Even with 
the complete activation of the cold walls, Hg’ will continue to persist at 
some level, as sources of previously accumulated mercury from previous 
thruster operation periods release Hg° which is, in turn, cryopumped on 
the baffles. Under maximum LN„ cooling of the shrouds and collector in 
the 5' X 11’ chanxber, total pressure (all gases) is 2 x 10 Torr, 
with principal fractions (assumed) of Hg®, H 2 O, N 2 and O 2 . 

The operation of the ion thruster will raise the measured 
chamber pressure some\»hat and also creates some neutral species that are 
not determined accurately by ion gauges (which, of necessity, must be 
somewhat removed from tho path of the thrust beam). A principal increase 
in ion gauge reading of chamber pressure due to thruster operation will 
be in Hg®, the bulk of which is probably at low energies ^^wall 
those portions of the chamber walls not refrigerated). In addition, the 
impact of energetic thrust ions on the collector and the shrouds creates 
the following gxoaaps of neutral particles: 

1) sputtered Hg® (previously cryopumped on locations impacted 
by thrust ions) , 

2) sputtered metal atoms from the collector and shrouds, 

3) sputtered H 2 O (previously cryopumped on locations impacted 
by thrust ions) and 

4) bouncing, neutralized, energetic Hg® (from Hg"*" impact on 
the collector with charge neutralization following, bttt 
not sticking). 

Of the species identified above, all may be expected to have energies 
significantly above wall temperatures during their initial flight 
(following sputtering and until an encounter with another wall). After 
the Initial flight, sputtered metal atoms will probably stick to any of 
the chamber walls and will remain accommodated there unless the area is 
subject to continued thrust ion bombardment. Groups 1, 3, and 4 above 
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may or may not stick to the first wall encountered after the initial 
flight. If the first wall encounter is an LN 2 cooled surface, sticking 
is expected. Accommodation and subsequent release at relevant vapor 
pressure rates would be expected for non-cooled portions of the chamber 
boundaries • 

The density of ”hot” (of the order of several eV for spattered 
particles) neutrals in testing chambers is non-trivial, considering here 
that sputtering ratios of Hg"*” thrust ions may range well above unity and 
that the streaming velocity of sputtered neutrals is significantly less 
than Hg"^ thrust ion velocity. In various regions of the chamber, then, 
weakly energetic neutral atom densities exceed the densities of Hg"*" 
thrust ions. As will be discussed in Section 3.4, these weakly energetic 
neutrals can react in specific forms to produce observable facility 
effect fluxes, principally in the backward moving ions. 

While weakly energetic neutrals are the principal cause of 
some facility effects, the major facility effect is that produced by 

thermal (wall temperature) mercury atoms acting in charge transfer 

“f* 

reactions with Hg thrust ions. Section 3.4 will discuss this reaction 
in further detail, and Section 5 will model expected fluxes of these 
charge transfer facility effect ions. 

3 . 3 Charged Particle Species 

Section 3.2 has noted charge transfer reactions between Hg^ 
and various neutral species in the testing chamber as a source of facility 
effects ions. In addition to these reaction produced ions, there will 
be small quantities of Hg"*" resulting from primary Hg"^ thrust ion bounce 
from the collector, without neutralization during the contact period. 

These backward moving (and probably weakly energetic) ions would not be 
distinguishable from weakly energetic (bouncing) Hg® which charge transfers 
against an Hg thrust ion to produce backward traveling ions. 

3.4 Facility Effects Reactions 

The measurements program reported here has been specifically 
directed to ion flux measurements, so that the facilities effects reactions 
of interest here are those which result in charged particles. As discussed 
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above, facilities effects ions can result from thrust ion bounce, without 
neutralization, upon impact with the collector. The expected 
bouncing ions is, however, small, and the major concern for facilities 
affects ions will be from charge transfer reactions between thruster 
ions and f acility.uienttals. In the charge transfer reaction. 


+ B" 


A“ + b'^ 


( 1 ) 


an initial ion. A^ and nautral, B”. exchanga aharga atata through electron 
tranefer. In assaaslng ralotlve magnltud.a to facllltlea effects, the 
created charge species will depend upon the combined densities of A and 
B- and the charge transfer cross section, The charge exchange cross 

sections attain largest values for "resonant" charge transfer. For 
example, the reaction 

Hg'*’ + Hg“ •> Hg® + Hg'*’ 


( 2 ) 


is resonant (AE - 0) and has a cross section of 'v 5 x 10 square 
centimeters, for relevant values of intraparticle collision energies. 

Of the several facility neutral particle species and thruster 
ion species, the two largest facility effects reactions both involve 
resonant charge transfer between Hg-^ and Hg®. In the first, the transfer 
occurs between Hg"^ thrust ions (Group I) and ambient chamber Hg®, large y 
at wall temperatures. This reaction perturbs measurements of "genuine" 
Group IV ions (Kg'", Hg® charge transfer downstream of the accelerator 
grid for ions and atoms leaving the thruster). The second major facility 
effect reaction involves sputtered (weakly energetic) Hg® charge transfer 
with Hg"^ thrust ions and causes backwards moving (0 ->180®) facility 
effect ions. A final facility effect of interest involves charge 
transfer between ambient chamber Hg® back diffusing into the ion 
acceleration space and producing a facility generated Group II ion which 
perturbs measurements of "genuine" Group II ion fluxes. 

While other forms of (a"^, B®) charge transfer occur in the test 
chamber (for example, sputtered metal atom, Hg"^ charge exchange), the 
relatively reduced magnitudes of non-resonant charge transfer cross 
sections indicates tha. these are not significant contributors to observed 

facility effects currents. 
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3.5 Assessment of Facility Effect Current Magnitudes^ 

Since facility effect ion currents will be measured with and 
will perturb genuine ion fluxes, it is desirable to assess the magnitudes 
of these spurious ions. In this program two approaches have been used, 
and each approach has limitations in the accuracy of assessment. The 
first approach, to be discussed in greater detail in Section 5.0, is 
analytical and involves calculations of ion generation between known 
thrust ion beams and modeled Hg® ambient densities. Uncertainties in 
this approach result from uncertainties in the ultimate deposition 
patterns of these calculated facility effects currents. A second, 
experimental, approach involves deliberate variation of ambient chamber 
pressure (for example, by controlling the temperatures of the LN 2 cooled 
shrouds) , observing ion currents as a function of measured chamber 
pressure and extrapolating observed ion fluxes to zero chamber pressure. 

The uncertainty in this approach results from possible variations in the 
partial pressures of the various chamber gases as total chamber pressure 
is varied. It should be emphasized, furthermore, that various gas species 
do not have identical ion gauge constants. An increase in ion gauge 
indicated chamber pressure by a factor of two, (for example, by allowing 
shroud temperature to rise) does not mean that ambient Hg pressure 
increased by a factor of two. Instead, a relatively higher fraction 
of the chamber pressure gain could have been obtained by increases in 
background water vapor than from Hg®, since cryopumping of H 2 O at the 
baffles loses effectiveness before a loss of cryopumping of Hg®. Thus, 
while chamber pressure variation measurements are interesting, there are 
possible inaccuracies, and, to reduce the perturbation of genuine ion 
flux measurements by facility generated effects currents, the best procedure 
is to have the maximum possible pumping of ambient Hg®. In the measurements 
to be discussed in the sections following, facility generated "noise" is, 
generally speaking, below the genuine ion "signal". As noted earlier in 
Section 3.1, however, for increasing physical separation and increasing 
angular divergence, spurious effects increase in magnitude relative 
to genuine currents, and, for certain energy, angular, and spatial 
regimes, become the dominant terms. 
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4.0 


LOW ENERGY ION MEASUREMENTS 


4.1 General Considerations 

Low energy ion measurements were obtained from retarding 
potential analyses of the probe signals of the 4*' Swinging and 
the Piggyback J_^. Low energy ion measurements were also obtained by 
the Weasel J_^ probes. While the probe data from the multi-collector Weasel 
probes is included in the Engine Operation Data volume, they will not 
be treated further in this report because of undetermined (and, perhaps, 
undeterminable) low energy ion trajectory refraction effects which may 
be present in the collection and analysis of these particles. 

The low energy ions to be treated in this section are 
primarily Group IV Hg"^ ions created by charge transfer reactions between 
Hg*^ thrust ions and Hg® atoms. If the Hg® atom in the reaction is in 
its initial passage through the test chamber (having emerged from the 
thruster and before the first encounter with a chamber boundary), the 
ion formed is a genuine Group IV and would be present for thruster 
operation in space. If the Hg® atom is an ambient chamber particle, 
the Group IV ion formed is a facility effect ion. 

The trajectories and velocities of Group IV ions (both genuine 

and facility generated) have been examined in References 3 and 4 and will 

be described further in Section 5.0 of this report. In brief, the energy 

of the ion immediately after the charge transfer is essentially the 

energy of the atom (^'^ kT where T is either the thruster wall temperature 

or chamber wall temperature). The ion is then acted on by the electric 

field structure in the plasma formed by thrust ions and neutralizing 

electrons. Potential variations in this thrust ion plasma are of the 

order of a few kT where k is Boltzmann’s constant and T is thrust 
e e 

beam electron temperature. For typical neutralization conditions, kT^ 
is approximately a few tenths of an electron volt, so that Group IV 
ions, moving in these electric fields, acquire energies which arc only 
of the order of electron volts and, hence, are easily altered in 
trajectory by electric fields in the sheath regions between the plasma 
beam and the surfaces of probes measuring these Ions. These refracting 
electric fields can severely perturb measurements of ion directionality 
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at the probe location. The deposition patterns of the Group IV ions 
not severely perturbed by the presence of probes, however, and, for this 
reason, data from the 4" probe (a wide acceptance angle, total ion 

diffusion current measuring probe) can be utilized. (Deposition pattern 
measurements are not affected because ion deposition points are determined 
by the integral of E along the path, and these Integrals and path lengths 
are substantially larger, in general, than for the Integral and path 
length of perturbation fields in regions surrounding the probe). 

One further consideration in the usability of J_^ probe signals 
for Group IV ion measurements is the extent of other, and competing, 
ion flux signals. For probe locations in increasingly dense portions 
of the thrust beam, the currents of Group I and Group II ions into a 
cup can be significantly larger than the low energy Group IV ions. 

Under these conditions the method of retarding potential analysis becomes 
increasingly subject to error from spurious effects arising from energetic 
ion currents. Bor this reason. Group IV measurements for large and 
increasing axial distance z are increasingly subject to error. A 
mitigating circumstance, however, is that these are not the regions of 
principal concern for Group IV deposition effects. 

^•2 Testing Chamber Ambient Pressure Effects 

Section 3.0 has discussed facilities effects and has noted the 
production of spurious Group IV ions by charge transfer between Hg"*" 
thrust ions and ambient Hg". Such effects can be examined by deliberate 
variations in the density of ambient Hg°, although (as noted in Section 
3.4) this procedure is imprecise because of possible variations in the 
relative partial pressures of ambient neutrals as overall chamber pressure 
is varied. 


To produce variations in ambient neutral density, the rate of 

liquid nitrogen feed to the upper shroud In the test chamber was varied. 

Under maximum LN2 cooling, chamber pressure readings from the ion gauge 
^ —6 

are in the range from 2 to 4 x 10 Torr. Reducing the LN^ feed to the 
upper shroud increases chamber pressure to values of 8 x 10”° Torr. 
These chamber pressure variations caused the variations in Group IV 
ion flux illustrated in Figures 9, 10 and 11. The probe in use for 
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Figure 9. Low Energy Hg”*" Charge Exchange Ion Flux In the 4" Probe 
as a Function of Axial Distance, z, for Two Testing Chamber Pressure 
Conditions for a 0.5 Ampere Thrust Beam. 
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Figure 10. Low Energy Hg*^ Charge Exchange Ion Flux in the 
4" J Probe as a Function of Axial Distance z, for Two Testing 
Changer Pressure Conditions for a 1.0 Ampere Thrust Beam. 







Figures 9 and 10 was the 4" used as a retarding potential analyzer 
to selectively record the low energy ion flux while the Piggyback J^, also 
used as an RPA, was employed for the data in Figure 11. 

From the data in Figures 9, 10, and 11 three conclusions may 
be drawn. The first of these is that facility effects are present in 
Group IV ion measurements. The second conclusion is that the effects 
of ambient Hg® in producing low energy Hg'*’ are not everywhere equal and 
that certain regions and certain directions of probe orientation are 
more affected than others. The final conclusion is that, within 
selected regions and probe orientations , and for reduction of chamber 
pressure into the 2 to 4 x 10 Torr range, the bulk of the observed 
Group IV ion flux is genuine, that is, has resulted from a charge 
transfer between an Hg"** thrust ion and an Hg® escaping from the 
thruster. Measurements of Group IV fluxes under these conditions, then, 
is a representative measurement of conditions that would occur for a 
thruster on a spacecraft. 

By extrapolating the data of Figures 9 and 10 to zero pressure, it 
may be seen that for regions near z 'v 0 for the 4" probe, most of the 
observed signal is due to facility presence. In the region from z 10 
centimeters to z ''' 25 centimeters, and for the best pumping conditions, 
the 4" probe signal is predominantly genuine Group IV. The evidence 
from the Piggyback J_^ is somewhat more complex. In order to enter this 
probe, an ion trajectory must be at divergence angles in excess of 

135® and, as will be seen later in the modeled ion trajectory calculations, 
it is difficult for genuine Group IV ions to attain these high backward 
angles. The bulk of the Piggyback J^, then, is probably a facility 
generated Group IV flux, for all z and even under the best of pumping 
conditions. It does appear, however, for z in the range of 10 to > 25 

centimeters that a fraction of the observed signal is the result of 
genuine Group IV production. It should be noted here, however, that cup 
currents to the Piggyback are substantially lower than the currents to 
the 4" J^. A part of this reduction may be attributed to reduced cup 
size and solid angle of acceptance for Piggyback compared to 4" .1^, 
but, since the total cup currents have relative values separated by 
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' 10 , at least a portion of this reduction must be due to absolute 
reduction in low energy ion flux magnitudes at high angles 135* to 180*) 
compared to their magnitudes in the 90* + 45* range. Thus, while some 
genuine Group IV flux may be present in the backward hemisphere, the 

magnitude is significantly lower than Group IV fluxes in the forward 
hemisphere. 

^•3 Characteristic Group IV Ion Flux Shane 

In examining the 4" determinations of Group IV ion flux in 
Figures 9 and 10 (and in the figures in the Engine Operation Data), a 
characteristic shape in the ion current pattern is apparent. The solid curve in 
Figure 12 illustrates this characteristic shape, identifies specific regions in 
the flow, and estimates the relative flux that would be obtained under 
space conditions where facility created ambient Hg* has been eliminated. 

In the regions of z < 5 cm (for the radial value of 28 cm 

of the 4" probe, at angles of divergence, 0 , greater than 80®) 
the probe current has an exponential behavior in z. In the range from 
z 5 cm to z % 20 cm, the exponential rise gradually lessens to a 
plateau at 20 cm. For z > 20 cm, the observed flux generally 
diminishes slowly before reaching a lower, almost steady, level. 

The dashed curve in Figure 12 is an estimate of the genuine 
level of Group IV Hg . In Region I, genuine Hg is dominated by facility 
effect ions. A similar situation arises in Region III, in which a second 
problem of high "noise" levels of Hg"^ thrust ions complicates the retarding 
potential analysis process out of which the low energy Group IV ions are 
Identified. Only in Region 11 is the signal-to-noise ratio acceptable. 

Probe location in Region II for determinations of Group IV Hg"*" 
as a diagnostic of thruster operation is appealing not only because of 
improved signal-to-noise ratios, but also because of the more gradual 
variance of ion flux as a function of z. (Location of diagnostic probes 
in high gradient locations is potentially troublesome in that only minor 
changes in plume "shape" can cause mf.jor flux changes at a given location. 

In the plateau portions of Region II, such effects may be expected to 
be considerably reduced). 
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Figure 12. Characteristic Hg'*’ Charge Exchange 
Probe and Estimated Genuine Flux as a Function 


Ion Signal in 4” 
of Axial Distance, z. 


Figure 12 's characteristic shape curves have been observed many 
times for the condition of both genuine and facility effect ions being 
measured. The estimated "true" flux (dashed curve) is subject to 
uncertainty. To fully remove uncertainty at z 0 by experimental 
measurements would require reduction in ambient chamber pressure to 
^ 10 Torr, requiring, thus, the largest available test facilities and 
a high level of cold wall activation. Eliminating facility effects in 
Region III also will require "v 10 ^ Torr chcunber pressure, noting, of 
course, that genuine Hg"*" signals fall off fcr increasing z in this region 
and that, ultimately, measurements for increasing z will become measure- 
ments of facility generated particles. 

It is also possible to reduce the uncertainty in the estimates of 
genuine flux by computing the expected genuine Group IV total production 
and comparing this value with along the cylinder on which the 

4" probe moves. Section 5 will discuss these computations further. 

For the present discussion it will be noted that ^J^dA from negative 
z values to z ^ 30 cm provides a total measured current which is in close 
general agreement with modeled total Group IV ion production. In 
Region II, then, and for chamber pressures remaining within the 4 x 10~^ 
Torr range, the signal-to-noise ratio probably remains above 1, so that 
a firm lower bound estimate of genuine Group IV Hg"^ in this region can 
be obtained by dividing the measured flux density by 2. Averaging 
measured and lower bound flux estimates leads to 

Group IV JJ % (.75 ± .25)JJ (3) 

genuine, Region II measured 

^•4 Slow Ion Behavior as a Function of Thrust Ion Current 

The 4" probe was utilized in a series of beam scans in which 
thrust ion current was varied. The Group IV ion current density was then 
examined as a function of z as a function of I^, using, as a first method 
of characterization the linear regression 

For most z the intercept term a(z) was generally small, leading to a 
conclusion that, under proper engine operation, the Group IV ion current 


27 


density scales, generally, ns thrust beam current 

. Table 1 lists the 

least squares fitted 

aj.y(z) and Gj,y(z) for z from 

-20 

cm to +40 cm at 

r 28 centimeters. 

Figure 13 illustrates values 

of ' 

4" cup current 

at various z values 

and as 1. is varied and the least 

X 

squares fitted 

linear regression of 

this data. 



Table 

1. Linear Regression of J_^ 

(Group IV) 

as a 

Function of Thrust Beam Current. 


z (cm) 

a^y(z) 


Gjv(z) 


(pA/cm^) 


(10"^cm"^) 

-20 

-.002 


.029 

-10 

-.009 


.095 

-5 

-.022 


.312 

0 

-.054 


.794 

5 

-.036 


1.462 

10 

.054 


2.184 

15 

.144 


2 ...52,7-. 

20 

.054 


2.383 

2j 

.686 


1.624 

30 

.360 


1.877 

35 

.902 


1.588 

40 

.415 


2.166 


The approximately linear behavior of with appears 

at first to violate reasonable assumptions of this dependence. Genuine Group IV 
ion production rate is clearly proportional to the product of the ion 
thrust current and the neutral atom current released by the thruster. 

If the thruster tended to operate at constant propellant utilization as 
I_^ varied, then 1^ (equivalent Hg“ release current) would die. proportional 
to I , and the I_^l^ product would be proportional to 1^. Some evidence 
of this dependence in Group IV ions had been previously observed with 
a 20 centimeter thruster (Reference 1). The 30 centimeter diameter I.eRC 
thruster, however, has been characterized as operating at Increasoc 
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propellant utilization for increasing thrust ion current, and Kaufman 
has observed mercury bombardment discharges in which remains essentially 
constant as varies, for optimized discharge operation. Accepting 
as fixed for varying would lead to the approximately linear behavior 

of Group IV ion current with l_j_, as given in Table I and Figure 13. 

A remaining question in these Group IV ion measurements is the 
dependence of facility generated Hg"*" charge exchange as is varied. 

For facility effect ions, total production rate will be proportional to 
the product of the beam current I and the ambient chamber density in 
Hg°. It would appear reasonable to assume that Hg“ ambient density is 
proportional to Hg^ beam current. For a purely cryopumped system, 
ultimate chamber pressure should be proportional to the rate of Hg"^ inlet. 

The observed chamber behavior, however, is that pressure does increase with 
increasing beam current, but the relationship is more of the form + 

k, I^ where k, ^ k-I. ^ (for I . _ 2 amperes). From this observed chamber 

1 4“r,t T,L 

behavior, the rate of facility effect Group IV ion production should be 

2 

proportional to k, I . ^ + k«I, . , and the relevant question then becomes 
the magnitudes of kj^ and k 2 to each other and to e(z) in Equation 4. 

An examination of the linear regressions of Table I in the range -20 < z < +5 
centimeters shows strong linearity. Since this range in z is expected 
to be significantly influenced by facility effect charge exchange ions, 
it would appear that these ions are also almost linearly dependent on 

l, . , and that the k, I term generally dominates. 

+ J t X T , C 

A final treatment cf slow ion behavior has examined the 
Piggyback signals as functions of z and I_^, again using linear 
regression. For a representation 

I_j_ (Piggyback) - Sp(z) + kp(z)I^^^ (5) 

the linear regression values given in Table 2 are obtained. In general, 
the probe signals were approximately linear in I, although more scatter 
in the data was experienced for this probe than in the 4" J_j_ data given 
before. 

The examination of the Piggyback data was not specifically 
to gain further insight into the behavior of genuine Group IV ions since, 
as discussed earlier, the principal signals at these high divergence angles 


will not be genuine but will, instead, result from the ambient chamber 
gas. The Piggyback data does, however, tend to confirm notions from the 
A" J_^ data that facility effect ions also scale linearly in 1 


Table 2, Linear Regression of Piggyback 
Signal as a Function of Thrust Beam Current. 


z (cm) 

HpCz) (nA) 

kp(z) (units of 10 

-13.5 

-6 

80 

-8.5 

-6 

101 

-3.5 

2 

110 

1.5 

18 

121 

6.5 

2 

203 

11.5 

-108 

546 

16.5 

22 

432 

21.5 

-160 

570 

26.5 

-30 

380 

31.5 

-36 

366 

36.5 

-28 

356 

41.5 

-120 

425 

46.5 

-60 

380 

Slow Ion Behavior 

as a Function of 

Screen and Accelerator 


Potential 

The Group IV ion flux was also examined as a function of net ion 
energy (screen potential varied for all other thruster potentials fixed) 
and as a function of accelerator voltage (for other thruster potentials 
fixed). Both of these variations lead to changes in the accel-decel ratio, 

R, where R = ^ l^gl^ where r is the deceJ-accel ratio used in 

later sections. 

For convenience here and for its use in later sections, Table 3A 
lists nominal thruster operation parameters at the various data points. In 
the data discu.ssed in this section, thruster operation will be at points 3 6 
and 18, 2 and 19, and 22, 23, 24. Table 3B lists actual engine operation 
conditions . 
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Table 3A, Engine Operation Data Points 
and Nominal Thruster Operational Parameters. 


CO 

0 

•H 

£ 

(0 


±.02 ±.l ±.01 ±.l ±.5 - ±.l 


1 1.1 2.0 0.5 12 37 Nom. 2.0 




25 1.1 2.0 0.5 8.6 43 Nom. 2.0 
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Table 3B. Thruster Operating Parameters During 
Beam Efflux Measurements. 


••m current 
^ propellant 
tllization 
^iciency 


it ion 
lergy 


:celerator 

Itage 


-Icharge 
; Itage 


utralizer 

■«per 

:rrent 


utralizer 
am rate 


to 

<u 


•H 

0 

(X. 

CO 

u 

CO 

p 

Screen 

voltage 

(kV) 

Beam 

current 

(A) 

: 

Accel 

voltage 

(kV) 

Emission 

current 

(A) 

Discharge 

voltage 

(V) 

Date 

Tank 

pressure 

(yTorr) 


Prop. 

utiliz. 

eff. 

1 

1.1 

1.94 

0.55 

12.0 

37 

10/29 

4.5 


102 

2 

3 

1.1 

1.90 

0.55 

10.0 

39 

10/30 

5.0 


90 

4 

1.1 

1.5 

0.55 

9.0 

37 

9/15,9/25 

3.8,5. 

7 

94,96 

5 

1.1 

1,46 

0.50,0.55 

7.5 

38 

10/30,11/21 

3.5,3. 

9 

74,78 

6 

1.1 

1.40 

0.55 

7.0 

39 

10/28 

4.8 


82 

7 

1.1 

1.0 

0.55 

6.0 

39 

9/10,9/11 

3.6,6. 

3 

84,86 

8 

1.1 

0.95 

0.55 

5.0 

37 

9/17 

3.4 


80 

9 

1.1 

0.91,0.94 

0.55 

4.5 

37,39 

9/17,11/21 

2.6,5. 

8 

62,79 

10 

1.1 

0.52 

0.55 
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Figures 14 and 15 Illustrate 4" current as a function of z 

for all ions, for soft ions (energy less than 25 eV) and for hard ions 

(energy greater than 95 eV) as thru.^ter screen potential is varied from 
1.5 kV to 0.7 kV, all other thruster voltages held fixed. Three features 

in the data given there are of interest. The first of these is that the 

Group IV ion flux is virtually unaffected for this screen voltage 
variation (and consequent decel-accel variation). This non-sensitivity 
of the Group IV flux will be shown (in Section 5) to be a reasonable 
consequence of expected trajectories for these slow ions. The second 
feature of interest is the observed increase in high energy, high angle 
of divergence ions as the decel-accel ratio is increased. This feature 
will be examined further in Section 6. The third feature of (general) 
interest is the relative magnitudes of hard and soft ion fluxes for z > 

10 centimeters. Sections 4.1 and 4.3 have discussed the problems of 
retarding potential analyses to determine low energy ion flux in the 
presence of large quantities of high energy flux (see also Figure 12) and 
Figures 14 and 15 illustrate the regions for which these low signal-to-noise 
conditions are obtained. 

Figures 16 and 17 illustrate 4" data, separated into low and 
high energy ions, for engine operation points 2 and 19, for which the 
accelerator grid potential is set at -550 volts and at -300 volts (a 
minimum value required to prevent electron backs treaming ) . The slow ion 
Group IV flux again demonstrates its lack of sensitivity to changes in 
thruster operation. A major change is experienced, however, in the high 
energy high angle-flux which has significant implications in terms of 
the use of Group IV as a thruster diagnostic (see also Section 8.0). In 
Figures 18, 19, 20, these data runs are repeated, except for a now lowered 
(1 ampere) thruster current. Again, slow ion flux is relatively invariant 
to changes in the decel-accel ratio, and, for minimum decel, there are 
significant diminutions in the high energy high angle flux. (Figure 18, 
Engine Operation Point 22). 

^•6 Slow Ion Behavior as a Function of Discha rge Potential 

The Group IV flux was also examined as a function of discharge 
chamber potential. Figures 21 and 22 illustrate the data for a 1 ampere 
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Figure 17. Total Ion Current, and Soft and Hard Ion Current 
Components in the A" J Probe as a Function of Axial Distance 
z, for Engine Operation Data Point 19 (V ■ -.3 kV). 
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Figure 18. Total Ion Current, and Soft and Hard Ion Current 
Components in the 4" J Probe as a Function of Axial Distance 
2 , for Engine Operation Data Point 22 (V - -.I kV). 
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Figure 19. Total Ion Current, and Soft and Hard Ion Current 
Components in the 4" Probe as a Function of Axial Distance 
z, for Engine Operation Data Point 23 (V - -.5 kV). 
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Figure 21. Total Ion Current and Soft and Hard Ion Current 
Components in the A” J. Probe as a Function of Axial Distance 
z, for Engine Operation Data Point 28 " ^3 V). 






thrust beam as discharge potential was varied from 43 to 34 volts (engine 
operation points 28 and 30 in Table 3) at constant discharge power. The 
results in Figures 21 and 22 demonstrate that the hard ion flux is 
essentially invariant to these thruster changes, and, that the low .energy 
ion flux is similarly unchanged. This behavior will be considered as 
logical in view of expected slow ion trajectory modeling (Section 5). 

4.7 Slow Ion Behavior as a Function of Neutr alizer Operation Conditio _ n 

The engine operation point data in Table 3 indicate variations in 
neutralizer conditions and the thruster plume was examined for two conditions 
of neutralizer keeper potential. In the first condition, neutralizer keeper 
was held at 13.6 volts. In the second condition the neutralizer tip heater 
power was reduced somewhat leading to a keeper potential of 17.0 volts. The 
4- J probe data for all ions and for Jiard and soft ion currents is given in 
Figure 23 for the two keeper potential conditions. As shown there, neither 
the hard nor the soft ion flux patterns are significantly affected by this 
variation in neutralizer operation condition. 

In separate experiments and using the Engine as a floating 

probe, evidence was obtained that, although the potential in the thrust 
beam did move upward as neutralizer heat rate was decreased, this potential 
increase did not lead to an observable increase in thrust beam neutralizing 
electron temperature. Since T^ did not vary significantly and since the 
thrust beam plasma density and density gradients were not affected by 
this variation in neutralizer operation, there is no expected basis for 
soft ion trajectory variation, and the data in Figure 23 confirms these 
expectations. This insensitivity in the Group IV flux to neutralizer 
operation is not expected to continue to be obtained, however, if severe 
cutback of neutralizer heat occurs leading to a modal change in neutralizer 
operation and significant increases in thrust beam neutralizing electron 
temperature. 

4.8 Slow Ion Behavior as a Function of Collector Surface ^ Materia l 

The discussion in Section 3 considered several forms of ambient 
chamber particles and it could be considered possible that the species 
and fluxes of these ambient gases would vary as the collector surface 


43 



material is varied. Although these possibilities were not pursued in 
detail, two conditions of collector surface were used in the beam 
measurements program. The first of these was a bare titanium collector 
(used here and in all other data points with the single exception of these 
experiments) and a low sputter yield graphite sprayed collector. 

The probe used in the ion flux measurements here was the 
Piggyback J^, since this probe looks directly at the collector surface 
and could be expected to respond more sensitively than any other probe 
to changes in the collector surface. Figure 24 illustrates Piggyback 
current for the two surface conditions. Although the ion gauge readings 
of chamber pressure were essentially the same for the two collector 
conditions, a significantly higher ion flux in the backward direction 
was obtained with the low sputter yield surface, compared to bare 
titanium. There is no immediate explanation for this observed behavior. 

As noted earlier, all other data runs were obtained with the bare 
titanium collector, whose use considerably simplified chamber operation 
during the thruster measurements (no required insertion and refurbishment 
of special collector surfaces). 

4.9 Slow Ion Behavior as a Function of Propellant Utilization 

Section 4.4 has discussed slow ion-behavior as a function of 

thrust ion current and has noted that propellant utilization varies as 

thrust ion current varies, with more efficient utilization of propellant as 

beam current increases. A statement of slow ion behavior as 1 varies, 

+ , t 

then, may Le that the Group IV production has I, ^ as an explicit variable, 

”T, t 

and propellant utilization as an implicit variable. 

It is also of interest to consider Group IV production with 
propellant utilization as an explicit variable. Such experiments are carried 
out by 1 aiding I fixed and varying propellant utilization by alterations 
in the bombardment discharge. In Section 5.2.3, total (measured) Group IV 
production is examined, and, within a given I, condition, various propellant 

X , t 

utilizations were maintained in the various engine runs. A qualitative 
observation is that Group IV production increases (generally) with diminutions 
in propellant utilization. 
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5.0 LOW ENERGY ION FLUX MODELING 

5 . 1 General Considerations 

The modeling of the low energy Group IV ion flux in the thruster 
plume can be carried out, for total charge exchange ion production, for 
the trajectories of specific ion;i ^following the charge transfer process) 
and for the flux deposition patterns of all ions. The analyses in this 
program have emphasized the first two areas above (total production 
and selected individual trajectories) but have not attempted to derive 
total flux deposition patterns. Discussion in subsequent sections will 
detail the reasons (increasing numerical complexity and diminishing 
levels of certainty in the analytical model for these calculations) that 
total flux deposition patterns have not been computed. The necessary 
steps leading to a total flux deposition computation will, however, be 
discussed. 

In addition to computations of genuine Group IV ion production 
and trajectory, the flux modeling has examined facility effect charge 
exchange ion production and trajectories for facility effect charge 
transfer involving either thermal (wall) neutrals or weakly energetic 
(sputtered) neutrals in the charge exchange. 

Some simplifications have been introduced in the model and in 
the calculations, of necessity. For example, neutral atom emission 
density from the thruster has been assumed to be uniform over the face 
of the thruster. Also Hg° emission from the neutralizer has been 
included (in terms of overall neutral release) but, rather than have this 
release emitted asymmetrically, the emission has been considered as a 
portion of the total release from the thruster discharge. 

A final point of emphasis here is in the complications caused 
by multiple thruster beams. These "cluster effects" will include an 
increased total production (since neutrals escaping from one beam without 
charge transfer may, in traversing the now adjoining ion beams, engage 
in such a transfer) and asymmetries in the total Group IV plume (at least 
within distances comparable to the thruster-to-thruster separation). 
Section 5.4 which discusses uncertainties in the modeling will also 
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exainine tha effects of increa;<ed Group IV production, including "pilo-up" 
of these slow ions in the interspace between beams and possible broadening 
of the divergence cones of these slow ions as a result of slow ion 
accumulation. 

^ Group IV I o n Production and Comparison to Observed 
Ion Flux 

^ 2 1 Calculated Genuine Group IV Production 

In this section the total charge exchange ion rate of creation 
will be determined between a total thrust ion current of I , in 
coulombs per second, and a thruster total neutral release, F^, in atoms 
per second. It will be assumed that all ions are Hg'*' and all neutrals 
are Hg®. The volume rate of genuine Group IV ion production, in 
coulombs per second per cubic centimeter is given by 


dn 


>cx 


dt 




a 

cx 


n 


ne 


( 6 ) 


where is ^hrust ion current density in amperes per square centimeter, 

*^cx ^^’8 ’ ^^8°) charge exchange cross section in square centimeters, 

density of Hg® in atoms per cubic centimeter for Hg® 
released by the ion engine. The total charge exchange creation rate is 
then given by the integral of Eq. (6) over all space. 

In setting the volume integral up for machine integration it 
has been useful to state the thrust ion and mercury atom density 
distributions in terms of normalized coordinate distances and certain 
functional plume shapes. Using a cylindrical coordinate system in which 
the z axis is the thrust beam axis and r = z = 0 is the center of the 
thruster face, the volume integral in Eq. (6) becomes 


4c. ■/ / 


( 7 ) 


where azimuthal symmetry has been assumed in both the thrust ion plume 
and in the thruster neutral plume. The coordinates r and z will be 
normalized to the thruster radius, b, by 
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"W" 


z = ub 

and 

r = nb 


( 8 ) 

(9) 


For this present experiment, the thruster radius, b, is 15 centimeters. 
This radius b must be distinguished from a second radius term, r^^, 
which will be used to denote the "core" radius of the thrust ion plume 
for the "parabolic core-exponential wing" thrust ion density model to 

be used here. In general, b ^ r^^. 

The thrust ion density distribution to be used has two principal 
regions. The first of these is a "parabolic core", and in this region 


J^,t(r,^) 


1 1 . 


( 10 ) 

core 




while in the exponential "wing" 

I. 




3,Cr^b + H”' 


2 




^obV / I r > r 


( 11 ) 

wing 

Ob 


The boundary to the core ragion at axial distance z is given by r r^^^ 

+ k z where is a term used to denote the rapidity of growth in r of 
thii core region for increasing z. Another term above, a^^^, is used to 
match the exponential drop-off density in the wing region. In the modeled 
calculations of the beam from the 30 cm thruster, the values r^^ = 10 cm, 
a = 5 cm and k = 0.2 have been used. This modeled plume has generally 
good agreement with the observed ion beam from this thruster and grid 
. «nd a . are also normalized to thruster radius using 

( 12 ) 

(13) 

The final computer model of given by 


c: u • w. 1 


Ob 


Ob 


r , = .b 

ob 


and 


a , = * b 

ob 
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(14) 


where p is a function which |.roduces the core and wing regions of Eqs. (10) 
and ( 11 ). 


The neutral plume density is given by 


n (r,z) = -—5 

ne 


TTb^V 


o , th 


|Y(a,n)| 


(15) 


where, yOx.ti) is a plume function for neutral release and is a 

thermal atom release velocity. In the neutral plume used for the computer, 
atom emission density is assumed to be uniform over the thruster face. 

The emission of a single source point into solid angle dh at divergence 

angle 0 is given by 


'^'^ne = k.cos^O 

dt 


( 16 ) 


In the calculations, three distribution forms (cos 6 , cos 0, and cos 0) 
were examined. The bulk of the computations were carried out for the cos 0 
release above, and the results reviewed here will be for that comparatively 
"broad" release pattern only. In principle, a more accurate fit to the 
neutral release might be possible by an expansion in terms of cos 0 , 
cos^O, and cos^O forms of release but, for present purposes, the simplified 
neutral plume appears to be adequate. Values of k^ in Eq. (16) are :ch 
that one integral of neutral release over the face of the thruster aud 
over 2 it solid angle of release directions yields the total neutral release, 

F , (in atoms per second), 
o 

When Eq. (7) is transformed, using Eqs. ( 8 ), (9), (14), and (15) 
tlie total charge exchange ion f<->rmation rate is given by 
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+CX 


'‘‘W ex \ 



ndnei j.r (• «^ i ^ ) r ( ' i ' ' ) 


(17) 


or 


•fex 
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wlujre c; denotes the total integral in r, and (* of the ion and neutral 
universal plume shapes. 


Values of have been computed as functions of the various 

parameters in Eq. (18). For convenience, the display of these calculations 
is in terms of propellant utilization rather than F^, since the usual 
description of thruster operation will be in terms of thrust beam current 
and propellant utilization efficiency. Figure 25 illustrates the expected 
total charge exchange ion production for ^ in the range from 0 to 2 
amperes and propellant utilization in the range from 70% to 90%. Other 
parameters used in the calculations given there are: b = 15 cm, r =10 cm. 

^ob “ ^ ^1 “ 'i" neutrals, o, = 5 x 10“^^ cm^, and 

'^o,th * ^ ^ cm/sec. 


For an ion thruster operating in the .5 to 2 ampere and 70% 
to 90% propellant utilization range, total charge exchange ion production 
may vary from v 1 to v 70 milliamperes. Section 5.2.3 will compare this 
calculated production rate with observed production rates derived by 
suitable integrals of the 4" slow ion signal. 

A final aspect of the computations to be considered here is 
the sensitivity of the integral in Eq. (17) to the form of the neutral 
distribution. The Integral G has been evaluated for cos 0. cos^O, and 
cos j release models for neutrals at fixed total neutral release with 
the parabolic core/exponential wing ion beam plume (r^j^ = 10 cm, a = 5 c 
^ 2 ^ “ *2) with the following results: 


G 


0.41 
0.46 
0.50 

As may be seen, the more narrowly di,stributed distributions lead to somewhat 
larger charge exchange ion production rates. The variances here, however, 
arc not considered significant in view of other and more basic limitation! 

in the model (specifically in the uniform neutral release flux as.sumptlon 
at the thruster face). 


f(^') 


me 
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2 . 
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5, J . Caicu Faci L i ( v i i r<.u:t^ i.V Frocluc tlou 

Calculation of tlie charge exchange ion production rate between 
till. I iirust l)oam and ambient Hg'“ is simplified if the comparatively 
reasonable assumption is made of a uniform ambient neutral density, 

For this coinUtion, the traver:,.J of each axial distance increment, dx, 
by n thrust ion beam I. creates 

+, t 


+CX 

— = n n 

dz +,t cx na 


(19) 


For an ambient Hg° density in thermal equilibrium with the 
upstream end of the testing chamber 20“C), an ion gauge reading 

of lO"^ Torr will correspond approximately to an ambient neutral density 
of ' 3 X 10^^ atoms per cubic centimeter. Using = 5 x 10 cm^ 
leads to 


dl 


+CX 


dz 


1.5 X lO"^ 1. 


( 20 ) 


■'■’V,:Torr Hg“ 

where, it should be emphasized, the only neutral density of significance 
is that of llg® (because of reduced charge exchange cross sections for 
other, non-resonant, transfer./;. It should also be emphasized, again, 
that the ion gauge responds to all chamber gases and that the partial 
pressure of mercury in these experiments is not known. 


From Eq. (20) it may be seen that Hg"^ charge exchange ions 
are produced at 150 uamperes per centimeter per ampere of Hg thrust 
ion current per ;;Torr of Hg® ambient density. For the region from 
z - 0 to z = 50 centimeters this would yield 7.5 milliamperes of Hg 

charge exchange per ampere of thrust ion current per cTorr of Hg^ 

ambient density. The ion gauge readings during the beam measurements 
avro in the range from 2.5 tc 5 Torr, some fraction of which Is not 
Hg' . for Kg'" densities ranging f 1 to 5 i.Torr, the expected production 
of iig/charge exchange current in the 0 ^ z 50 cm interval would range 
from 7 to 21 m j 1 1 [amperes lor each ampere of Hg^ thrust ion current. 


The facility effect. Hg charge exchange ions may be compared 
to ih»i calculated genuine tiroup IV' in Figure 25, and may l>e st'en to he 




PROPELLANT UTILIZATION EFFICIENCY 

(PARABOLIC CORE - EXPONENTIAL WING IONS; 
r^b ^ 10 cm; - 5 cm; =.2) 

Figure 25. Computed Total Hg"*" Charge Exchange Current Formation as a 
Function of Thrust Beam Current and Propellant Utilization. Ton Beam 
is Parabolic Co re/ Exponential Wing and Neutral Emission is Uniform 
Over Thruster Face and Cos h Angular Distribution. 
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comparable in magnitude. This would tend to indicate, in turn, a 
comparatively low signal“tO“noise ratio in the determination of the 
genuine charge exchange ions. There is, however, some expected relief 
through concentration on measurements of Group IV ions in specific axial 
intervals. Section 4.3 in discussing the characteristic shape of the 
charge exchange current as seen by the 4" identified three regions 
(see also Figure 12) with Regions I and III having a greater dependence 
on ambient effects while Region II is comparatively well detennined by 
the presence of genuine charge exchange ions. 

As another aid to assessing comparative magnitudes of facility 
effect and genuine charge exchange currents, the fraction of all genuine 
Group IV ions created in the axial distance interval from 0 to z has been 
calculated. Figure 26 illustrates the value of the volume integral in 
Eq. (17) from 0 to a (z = cxb, where b is thruster radius) compared to 
the volume integral from 0 to •». It may be seen that 60% of the genuine 
charge exchange production occurs within the interval from z ® 0 to 
z = b (= 15 centimeters) and that some 78% of the genuine production occurs 
in the interval from z = 0 to z * 2b (30 cm) . By concentrating the 
measurements on axial regions near the thruster, good signal to noise 
conditions can be obtained. For measurements away from these regions 
(and as discussed in Section 3.1) the effects of facility presence are 

Another means of viewing the importance of various regions in 
the chamber is the calculation of the boundary along which the density 
of ambient Hg“ is equal to the density of Hg“ escaping from the thruster. 
Figure 27 illustrates the boundary in z/b and r/b for which n^^ » n^^ 
for an ambient density of 3 x 10^° Hg® per cubic centimeter and for an 
equivalent Hg* release of ^ 280 milliamperes (a propellant utilization 
of 0.78 at 1 ampere of Hg^ thrust current, and 0.88 at 2 amperes of 
thrust ion current) of "cos ' " neutrals at 500®K. Along the Indicated 
boundary the production rates of genuine charge exchange and facility 
effect charge exchange are equal, and, within the boundary, genuine 
effects predominate. Increasingly so for regions nearer and nearer the 
thruster face. 
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A final aspect in the comparison of genuine to facility effect 
ions to be noted here is that volume rates of production and deposition 
patterns of these ions are linked by comparatively complicated ion 
trajectory factors which are not the same for the two forms of production. 

These trajectories will be discussed further in Section 5.3. The use of 
the variances in trajectories can improve signal-to-nolse ratios in the 
measurements by selecting certain positions for measurement, as in 
Region II of the 4" probe movement. 

5.2.3 Comparison of Observed Group IV Ions to Calculated Production 
Rates 

The charge exchange ions produced by (Hg , Hg®) charge transfer 

will move both radially and axially. If their movement is predominantly 

radial (and Section 5.3 will demonstrate this predominant direction of 

motion), the ions should be capable of a straightforward measurement by 

the 4" probe. The integral of the current density of Group IV ions 

seen by the 4" probe over an appropriate range in z (and assuming 

azimuthal symmetry in the Group IV flux) would then provide a measurement 

of both genuine and facility effect Group IV ions. Table 4 contains 

values of ^ 2irr dzJ^ for a series of measurement conditions. The 
0 P 

measured production given there may be compared to calculated production 
(Figure 25) and estimated facility effect production (Eq. 20). Consider, 
for example, the data obtained at Engine Operation Point 12, an ion thrust 
beam of 0.5 ampere. From Table 4, the value of J* 2TtrpdzJ_^ from 0 to 50 cm. 
in z is 14.1 milliamperes. From Eq. (20) and for an assumed level of 
2 uTorr of Hg® ambient pressure, the expected facility production of 
Hg**" charge exchange is 7.5 milliamperes in the range from 0 to 50 cm 
in z. This would tend to indicate a genuine production of Hg of 

6.6 milliamperes. From the relationships in Figure 25 it may be seen that 
an 0.43 ampere thrust beam from a thruster operating at 52% propellant utiliz- 
ation efficiency could produce this indicated level of genuine charge exchange 
ions. There’ is, thus, a qualitative agreement between observed charge 
exchange production (from all sources) and expected production from 
genuine causes and from facility effects. 

While the nbove comparisons are comforting in terms of signal- 
to-noise ratios for ions over the total 0 < z < 50 cm interval, it should 
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Table 4A. Engine Operation Data Point with Thrust 
Ion Current, Propellant Utilization, and Calculated 
Total Genuine Charge Exchange Ion Current from Model. 
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Also given above is calculated total genuine plus facility measured ion 
current from 4 signal and integral 0 to 50 cm in z. Measured total 
production and calculated genuine production infers the facility effect 
generation and infers a facility partial pressure of Hg°. For internal 
consistency, the Hg — pressure may be compared to ion gauge reading of 
total pressure. Internal consistency good on lower range (.5A to lA) 
of thrust ion current. 
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Table 4B. Engine Operation Data Point with Tfirust 
Ion Current Propellant Utilization, and Observed 
Soft Ion Current in 4" Cup at z = 20 cm. 


Engine Operation 
Data Point 

Facility Total 
Pressure (viTorr) 

Maximum Permissible 
Soft Ion Energy in 
4" Cup (eV) 

— 

Thrust Ion Current, 

1 

Bombardment Discharge 
Emission Current (A) 

Propellant 

Utilization 

Soft Ion Current 
(in uA) in 
4" Cup at 

z = 20 cm 

1 

4.5 

25 

1.94 

12.0 

>1 

f 

200 i 

2 

5.0 

25 

1.90 

10.0 

.90 

250 

3 







4 

5.7 

20 

1.50 

9.0 

.94 

180 

5 

3.9 

25 

1.46 

7.5 

.78 

270 

6 

4.8 

25 

1.40 

7.0 

.82 

260 

7 

6.3 

9 

1.00 

6.0 

.86 

150 

8 







9 

5.8 

25 

0.94 

4.5 

.79 

240 

10 

3.0 

10 

0.52 

3.0 

.80 

46 

11 

4.2 

10 

0.50 

2.5 

. 66 

78 

12 

3.7 

25 

0.45 

2.0 

.54 

86 


For fixed thrust ion current and varying propellant utilization, genuine 
charge exchange ion signal should vary as n (l~n) where n is propellant 
utilization fraction. Observed results indicate rising production rates ^ 
for diminishing propellant utilization. Three factors restrict these 
conclusions to a-qualitative statement: 1) the range of propellant 

utilization variation, within a set I.f,t* generally not large, 

2) inaccuracies exist in propellant utilization measurements (reduced 
throughput measurement times and (possible) inventory fluctuations), 
and 3) facility effect variations also occur from one engine run to 
another. It should be noted, also, that engine operation throughout the 
experiments maintained generally good utilization (for a given 
that deliberately "spoiled" (greatly lowered) vitilization did not occur. 
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be emphasized again, that certain regions are expected to have a 
predominance in facility effect currents. Figure 12 has given these 
regions, and Region III is an expected facility dominated zone. The 
integral of fluxes in Region III for the Data Point 12 example above 
also shows a comparatively large fraction of the integral (50%) occurs 
in the range 30 < z < 50 cm. In this region, the observed signal is 
certainly facility dominated with a low signal-to-noise ratio. In the 
region from 10 to 25 cm in z, however, the observed currents are 
probably genuine for the most part , confirming the earlier estimate that 
the use of observed J_^_ in this range is probably an over estimate of 
the genuine charge exchange flux by only about 30%. 

^•3 Charge Exchange Ion Trajectories 

5*3.1 Genuine Group I V Io n Tra jjectorles 


The procedure for the calculation of charge exchange ion 
trajectories in the plasma thrust beam internal electric field structure 
was developed in Reference 3 and Reference 4. In this calculation the 
potential in the thrust beam is assigned the form 


V(r,z) 



P (r,z) 

Pq(0,0) 


( 21 ) 


where T^ is thrust beam neutralizing electron temperature, k is Boltzman's 
constant, and e is electron charge. For kT^ expressed in electron volts, 
e is assigned the value of unity, and V is in volts. The plasma density 
is p (ions/cm ) and is thrust beam density at r = z = 0, which is the 

maximum density in the parabolic core/exponential wing density model 
given in Eq. (10) and (11) and the potential at r = z = 0 is assigned the 
value 0. This electrostatic" barometric equation has been determined 
in Reference 10 for plasma thrust beams and found to be generally 
adequate in radial scans of the potential. There are questions on the 
adequacy of Eq. (21) for variations in the position of a field point 
in the axial direction and these will be discussed further in Section 5.4. 


From Eq. (21) and the density model in Eq. (10) and (11), 
variations in potential for movement Ar at (r,z) and Az at (r,z) have 
been calculated by a computer program. The value of E (r,z) and 
E^(r,z) are, thus, knoivn at all (r,z). 
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The remaining element in the calculation of charge exchange 
ion trajectory is a specification of the initial charge exchange speed 
and direction. For genuine charge exchange ion production, initial ion 
velocity is the atom velocity prior to charge exchange and is determined 
by the thruster "wall" temperature, the point of exit of an atom from 
the thruster and the point (r,z) at which charge exchange occurs. These 
three quantities are programmed into the calculation which derives the 
initial f and z of the charge exchange ion. 

The motion of the charge exchange ion is determined by a 
"marching" method in which particle acceleration from and E^, and 
known r and z at the point (r,z) is used to calculate the (r*,z*) and 
(f',z') of the ion after a passage of 6t in time. The electric fields 
at (r',z') are then calculated and the process iterated, with an 
accompanying trajectory printout. 

Figures 28 through 36 illustrate a series of trajectory plots 
of these charge exchange ions. Assumptions of thruster wall temperature 
(500°K) and thrust beam neutralizing electron temperature (5000®K) are 
stated on the figures. A particular simplifying limitation in these 
calculated trajectories is that the atom exit point and the charge 
exchange point occur in the same plane as the thrust beam axis, so that 
there is no ^ term and the calculations are only in two dimensions, r and 
z. The neglect of angular momentum is deliberate, and Section 5.4 will 
discuss the limitations which this restricted case places on total flux 
deposition patterns. 

The results of the trajectory calculations are not easily 
characterized. Trajectory can be seen to depend sensitively on the 
point of the charge transfer and on the direction of motion of the atom 
at the instant of transfer. Ions created in one side of a beam can cross 
the axis, in some cases, and emerge on the opposite side. Undoubtedly 
there are even more complications when angular momentum is admitted and 
the particle then has motion in r, z, and {/. 

In spite of the complexity, however, some general observations 
may be made. The first of these is that genuine charge exchange ions do 
not emerge in backwards directions for angles significantly over 90“. 
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Figure 28. Computed Hg^ Charge Exchange Ion Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point. 



CONDITION I ION BEAM i 

= ,67b, Oob = .33b, = .2) 

T A =10 

e w 

ENGINE NEUTRALS ( • = SOURCE POINJ) , 



3'- 


0 1 2 3 4 5 

(*/b) 

Figure 29. Computed Hg^ Charge Fjcchange Ion Trajectories for Knglnc 
Released Neutrals at Indicated Source Point and Charge Transfer Point 
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Figure 30. Computed Hg^ Charge Exchange Ion Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point 
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I i^ure 34.^ Computed l(^,"^ Charge Exchange ton Tra jector 1 cs for Engine 
Released Neutrals at Indicated Source Point and Ctiarce Transfer Point 
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" Lgure 36. Computed Hg^ Charge Exchange ton Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point 



Some ions do acquire the backwards hemisphere divergence cone, but the 
g. lat majority are directed '.utward at smaller divergence angles. A 
second observation is that ions formed at z generally would be perceived 
by a probe like the 4" probe at larger z values. This can be used 
to explain the general rise in observed prole signal in Region I (see 
Figure 12) to its plateau levels in Region II. 


One additional series of trajectory calculations were carried 


out for a modified thrust beam shape. In this new beam shape, the 
divergence constant, k^, in Eqs. (10) and (11) was reduced from its 
value of .2 in the examples above, to .1. This more narrowly diverging 
beam was expected to possess larger radial fields (compared to the 
axial fields) than for the more broadly diverging ion beam (k^^ - .2). 
Increased divergence in the electric field angles in the k^^ «^.l beam 
should, in turn, lead to increased charge exchange ion divergence. This 
"counter" motion between the thrust ions and charge exchange ion plumes 
(increasing divergence in one leads to diminished divergence in the 


other, and vice versa) is expected from general considerations and, from 
the results in Figures 37 - 41, is confirmed by the calculations. The 
observed effect is not a major reshaping of the charge exchange ion 
plume, which is a significant result if alteration of the thruster grid set 
or electrode voltages leads to an alteration in the thrust ion divergence cone 


^ ^ ^ Fac ility Effect Charge Exchange Ion Trajectories 
5. 3. 2.1 Thermal Atom/ Charge Exchange Ions 


The calculations of charge exchange ion trajectories have been 
carried out for two forms of ambient facility atoms. In this section, 
thermal atoms (T^ 'v 250»K) will be examined. The value of T^ has 

remained at 5000°K so that the temperature ratio, T /T , is 20 for these 
calculations. ® ” 


The principal difference for thermal atom/charge exchange ion 
trajectories from the engine atom case is that, for ambient atoms, all 
atom directions, prior to the electron transfer, are possible. A backward 
moving atom, thus, could continue in the backward direction after charge 
transfer, and would appear as a possible source of the previously observed 
facility effect ions (see Region I of Figure 12) near z - 0 for the 
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Figure 37. Computed Hg^ Charge Exchange Ion Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point 
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Figure 39. Computed Hg"*^ Charge Exchange Ion Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point. 
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Figure 40. Computed Hg^ Charge Exchange Ion Trajectories for Engine 
Released Neutrals at Indicated Source Point and Charge Transfer Point 




Figure 41 Computed Hg'^ Charge Exchange Ion Trajectorlee tor Poo. 
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4" probe. 

While these back diffusing ions are possible in principle, 
the calculations have not revealed strong tendencies for such motion 
for the thermal atom facility effects. Figures 42 and 43 illustrate 
trajectories for charge transfers occurring at four selected ^joints In 
the plasma thrust beam. Although some backward ion trajectories are 
observed, the flow tendencies are more accurately described as being 
radially outward. 

The results of Figures 42 and 43 do not appear to satisfy 
the observed shape of the charge exchange plume and two principal 
reasons may be advanced for this failure. A first reason is to note 
that the backward ions may be the result of charge transfers with more 
energetic ^atoms (Section 5. 3. 2. 2). A second possible explanation is 
that the E fields in these "wing" regions are not accurately described 
by Eq. (21) and the gradient calculation process. As Section 5.4 will 
note, the presence of significant levels of charge exchange ions in the 
dilute outer edges of the thrust ions invalidates plasma density gradient 
calculations based solely on thrust ions. The E fields in these wing 
regions could^be perturbed, with the perturbation acting to force the 
direction of E into the backward direction. These factors will be 
discussed later under charge exchange ion "pile-up" effects. 

Weakly Energetic Atom/ Charge Exchange Ions 

Section 3.2 has noted that the ambient atoms also include a 
weakly energetic component group from the sputtering actions of Hg'^ 
thrust ions on cryocollectors whose surfaces have resident Hg“. The 
energies of these sputtered atoms could be in the range from fractions 
of an electron volt up to electron volts, and, because of these larger 
atom energies, the trajectories of the ions following charge transfer 
can he continuations of these backward motions. Figure 44 illustrates 
calculated charge exchange ion trajectories for atoms at 5000“K (which 
IS also Tg, so that T^/T^ *1). The calculations clearly Illustrate that 
this increase in atom energy results in more widely dispersed deposition 
points for those atoms entering into charge transfer reactions. These 
trajectories may help to explain portions of both the 4" signal 
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Figure 43. Computed Hg"*" Charge Exchange Ion Trajectories for Ambient 
Chamber Neutrals at Indicated Charge Transfer Point and Initial Atom 
Direction. 





(near z = 0 for that probe) and the Piggyback signal (at all z) . For 
concinued increases in atom energy, the resulting charge exchange ion 
trajectories will become increasingly more straight line continuations 
of tlie earlier atom direction of motion. 

5.3.3 Comparisons of Observed Group IV Plume Shape to Calcul a_ted 

Trajectories 

5 . 3 . 3 . 1 Total Ion Flux as a Function of Axial Distance at Fixed _r _=_Tp 

Figures 28 through 44 illustrate calculated charge exchange 
ion trajectories for both genuine charge transfer (with atoms from the 
thruster) and facility effect transfer (with ambient chamber atoms). 

Figure 12 illustrates the characteristic plume signal from the 4" 
with an approximated "genuine" signal and Figures 14 through 23 lll'j.strate 
4" J data for various engine operation conditions. An important question 
is the internal consistency between these measured and calculated quantities. 

The calculations have shown that charge transfer to an Hg® 
escaping from the engine occurs, for the most part, within 2b 
(b = thruster radius) from the engine face, and that the ions so formed 
will generally encounter the cylinder which the 4" is 

located) at z values larger than the z at which charge transfer occurred. 
These calculated trajectory features explain the sharp rise in 4" 
probe signal for z > 0 and increasing. In many of the 4" J_^_ probe signals 
(Figures 14 - 23), this current reaches maximum levels near z 15 to 20 cm 
and has a perceptible (though small) decline for further increases in 
z to - 30 cm. This falloff in probe signal can be attributed to a falloff 
in genuine Hg'*’ formation for increasing z moving away from the thruster. 

If the ambient chamber density were zero, it is expected that the observed 
falloff would be more evident. The facility effect ions are, however, 
numerous and their presence leads to a continued high level of 4 
probe signal in the range above z 20 cm. An improvement in facility 
pressure to the range below 1 ,.Torr would, it is felt, clearly reveal the 
dropoff of ion signals in these downstream regions. 

The Piggyback piv'be signals (Figures 11 and 24) are not 
explained by calculated ion trajectories of genuine Hg Group IV or of 
Group IV from thermal Hg® in the test facility. The calculations of 
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trajectories for llg"*" formed from weakly energetic C'- eV energies) 
sputtered atoms (see Figure 43) can serve to explain some of the Piggyback 
signal. Since the currents here are small, however, (fractions of a 
microampere), other causes must be examined. These other signal sources 
include low energy Mg'*" trajectory refraction into high angles of 
divergence by electric fields in the sheath regions betv/een probe surfaces 
anc the charge exchange ion plasma, and a general broadening of the cone 
of divergence angles of charge exchange ions from "pile-up” effects 
(Section 5.4.1) . 

5 . 3 . 3 . 2 Angular Dispersion Pattern of Group IV Ion Flux as a Function 

of Axial Distance at Fixed r = r 

p 

The angular dispersion patterns of Group IV ions are determined 
by the Weasel I and Weasel II probes, which are mull.ico Hector retarding 
potential analyzers, moving at fixed r (^ 30 cm) in the axial range from 
-10 to +30 cm. These probes are illustrated in Figure 8 in Section 2. 

Data from the Weasel I and Weasel II probes is given in the 
Engine Operation Data volume, and will be discussed in this section in 
qualitative terms. One reason for this reduced emphasis on Weasel probe 
data has been previously noted, i.e. a concern that the ion flux densities 
and comparatively long flight paths within the analyzer, combined with 
the small value of Group IV ion energy, may result in ion trajectory 
refraction. A second reason for de-emphasis of the Weasel data are 
signal-to-noise effects in the accumulation of the retarding potential 
analysis data. As-noted earlier in discussions of RPA in the presence 
of strong hard ion current backgrounds, the presence of secondary electron 
emission from the collector by energetic ion impact can create the 
appearance of a "negative" population in the soft ion spectrum. For 
example, in performing a retarding potential analysis, the retarding 
potential grid is moved upward from 0 volts to a small positive voltage 
to prevent the passage of soft ions. Under normal conditions this leads 
to a diminution of collector current as the soft ion flux is no longer 
present. If secondary electrons are present, however, the positive; 
motion in potential of the retarding potential grid can cause tliese 
electrons (which were previously suppressed) to now leave the collector 
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surface to be collected at th.' i^P grid. This /reln.ase of eJ ec trua.-^ by 
the collector has the appearance of an increasing positive current 
signal at the collector, leading to the above mentioned "negative" 
population of soft ion states. 

Spurious signals in the soft ion detection, may be expected 
whenever the ratio of hard ion flux to soft ion flux is about equal to the 
inverse of the coefficient of secondary electron emission, Since 6 
for clean surfaces under energetic Hg"^ impact is generally small, (-y .1), compared 
to unity, a condition where j.) > 1/A requires a significant 

level of hard ion flux. From the 4" probe data, such regions are 
encountered at large z. For the Weasel probe this condition is 
encountered over a wider range in z. The reason for the difference 
in ranges of sensitivity is that the 4" probe has a large solid angle 
of acceptance, which tends to increase soft ion collection (broadly 
distributed in direction) over hard ion collection (most of which 
originates at the thruster face and emerges in a comparatively narrow 
divergence cone of directions at any one location). The Weasel probes, 
however, have a small solid angle of acceptance and soft ion current 
collected at any one collector is reduced considerably from the levels 
seen in the 4" J_j_. The net result is a significant impact on the 
quality of the Weasel probe data compared to that of 4" 

While the Weasel probe data is affected by these considerations, 
it has been included in the Operation Data Volume and will be reviewed 
here briefly. Near z = 0, the soft ion flux is chiefly present at 
0 = 90° and in the two backward channels (0 = 105° and 0 = 120°). From 
other discussion and calculations, a major part of the signal at z = 0 
is facility generated and might be expected to be directed into the 
backward hemisphere. For z in the 10 to 20 cm region (Region II of 
Figure 12), charge exchange ions appear in all channels, with principal 
concentration near 90°. The flow, thus, is generally a radial flow. 

While some improvements may be made in the Weasel probes, the 
procedure recommended for future measurements of the dispersion of tills 
soft ion "cloud" is to allow the flow to proceed without perturbation 
to separation distances which are large compared to the general scale 
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of tlio charge exdmnge production region. Such condJtlons can be 

. , Ifln< lpal, In ,i Urge u-stlng d.m,bnr »ltii extremely low 

mblen. mercury .eneity. II. e use ol Urge testing e„nm,.er wit, .nut cry 
w ucKgrounl llg- g tty U not sutlleient. since, n, nlreugy noted 
-lUty genernted llg u„s cun possess u seele sUe comparable to tbe 
Member length and there is no possibility of separation from these 
spurious currents or of probe placement at distances Urge compared to 
ion production volume scale size, 

Umltau^'s^o^ 

of Thrust Beam In ternal Potential Mn Ssi 

ass a I “'“bse exchange Ion trajectories have 

E<l' (21), the plasma beam Internal potentia:. 
rmu ation, which provides, via the appropriate derivatives, the 

and unctions which govern charge exchange Ion motion. There are 
wo presently recognlred limitations In this potential formulation. 

of -h h l‘>"itbtion in Eq, (21) appears to be present IrrespectJ 

-be Charge exchange Ion build-up and Involves the axial potential 

gu rent. Ue form of Eq. (21). which Is an electrostatic modification 
he baromutrtc equation, assumes that electrons of temperature T 
woul populate the plasma according to the hnown Ion density dlstrlLtlon, 
|^+,i(r,r). but without any net diffusion, while the net diffusion of 
elec tons in the radl^ direction Is essentially aero (negleetlng here 
comparattve.y small thrust Ion radial spreading), electrons must 
ompany the axial motion of the ions (v^:. . ) which calls for a 

Of E derivatives 

<1. (21). This uimlnution of axial electric field should result In 

argsr divergence angles for the total Internal electric field, 

for h " divergence angles 

for charge exchange ion crajectorie.s. 

the thrust b attempts were m.ide to determine 

. t her. internal potential using the Engine .1. probe (both using 

bb collector) as a floating prohe. w,.He these floating potltla, 

bsuremunts did reveal muasurahU r.adUl electric fluids, potuutlal 

variations In the axial direction were sufficiently small to avoid 
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detection, even for comparatively extei sive , (0<z<50 cm), axial motion of the 
floating probes. These measurements would appear to confirm the notion of axial 
potential gradients smaller than predicted from Eq. (21) derivatives. 

The radial electric field, however, is probably described accurately by 

the computation program, Eq. (21), and the known thrust ion density 
distribution. 

While the conclusions above are valid in the limit of the 
thrust ion density as the predominant ion density, another limitation 
(the second of the two noted earlier) arises from the growth of the 
charge exchange ion density to that point that the plasma ion density 
term in Eq. (21) is no longer given from the known thrust ion distribution. 
Charge exchange ion density growth, or "pile-up", would appear, at first, 
unlikely, since the currents of these particles are small compared to the 
thrust ion current. It should be noted, however, that the density of an 
ion species is describable as a current density divided by the ion flow 
velocity and that charge exchange ion flow velocities are approximately 
two orders of magnitude less than thrust ion velocities. Charge exchange 
ion densities can, thus, become comparable to thrust ion densities, with 
resulting perturbation to Eq.,(21). _ 

An estimate of conditions under which charge exchange ion density 
can become comparable to thrust ion density can be derived by considering 
a cylindrical column of radius r of thrust ions of constant current density J 
moving at v_^^^ through a neutral gas density of n^^. The volume rate 
of charge exchange ion production 


~**cx ^ coulombs . . 
dt ^sec cm3 ^ “ “^+,t ®cx% 

given in Eq. (6) predicts a total charge exchange production in a 
cylindrical volume of length, d£, and radius r of 


( 6 ) 


dN 


+CX 


dt 


irr d:l(J 


+,t 


a n ) 
cx n 


( 22 ) 


These ions, diffusing radially, will remain within the reference volume 
for a dwell time roughly given by 


1 
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+CX 


(23) 


leading to charge exchange ion density given approximately by 


D 


, a n \ 
+ , t cx n I 


(24) 


This charge exchange ion density will have a value relative to p 

"T, t 

(the thrust ion density) of 


‘^+cx^‘^+, t 


(0 n r) 
V . cx n ' 

+CX 


(25) 


2 

As noted earlier, the velocity ratio v, ,,/v is 10 , and 

p..„/p. ^ becomes of order unity for a n r of order lO”"^. Using 
C3C 1 ' ) c 13 2 cx n 

0 ■ 5 X 10 cm and r (for example) ® 10 cm leads to a requirement 

on neutral density of 

n^ < 2 X 10^^ atoms/cm^ (26) 

for the given choice of radius. The density given above 

in Eq. (26) corresponds to a chamber pressure of Hg” of 6 uTorr. Chamber 
pressure conditions of several pTorr, thus, can lead to charge exchange 
production rates sufficient to produce pile-up effects. 

The formulation given above is clearly an approximation and a 
complete and rigorous treatment of pile-up is beyond the scope of this 
program, since charge exchange ion presence alters the electric fields 
causing charge exchange ion diffusion, requiring, ultimately, a 
completely self-consistent iterative calculation of the formation and 
dispersion of these low energy particles. It is possible, however, to 
perform a slightly more detailed calculation of these effects. Figure A5 
illustrates plasma beam density using Eq. (21) and a parabolic core/ 
exponential wing ion density model, allowing charge exchange formation 
according to Eq. (6), and radial diffusion only. In the approximation 
there charge exchange ions move at kT^, where T^ is the wall temperature 
of the atoms, until that radial position in which 0 ^ “ Pj. Inside 

this radius the plasma potential is set "flat" (no electric field). 
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Outside of the radius, the charge exchange ions are assigned velocities 
which are self-consistent with Eq. (21), and the calculated ^ 

The curves in Figure 45 illustrate increasing perturbation? 
to plasma internal density and potential structure for increasing levels 
of Hg® pressure. Clearly it is desirable to maintain as low as possible 
an ambient density. 

Two final aspects of the calculations here are the multi- 
dimensionality of the effects, and the possibility that pile-up of these 
charge exchange ions can cause the "refraction" of charge exchange ion 
trajectories into the backwards hemisphere (as, for example, appears at 
the Piggyback probe) . The first aspect noted above is that the effects 
of charge exchange ion presence become more marked for diminishing 
thrust ion density. Thrust ion density diminishes in distant axial 
locations and in the exponential wing regions. Thus, perturbation of 
the internal potential of the plasma thrust beam is expected for large 
z, even on axis, and for large r, even for small z. This clearly affects 
many of the regions treated in the earlier trajectory calculations. 

Effects here are obviously complicated and no definite estimate can be 
given to the variation of charge exchange ion trajectories. One possibility, 
however, is that charge exchange ion build-up, particularly at large z, 
causes sufficient elevation of plasma potential at those locations to 
cause charge exchange ions formed at upstream locations to redirect their 
flow to higher divergence angles in order to find less thoroughly 
populated escape routes. 

5.4.2 Particle Coordinate Description Limitations 

The charge exchange ion trajectory calculations for Figures 28 
through 43 were carried out in two dimensions (r and z) only, neglecting 
azimuthal velocity, '|i, because of the complications which the inclusion 
of this coordinate would nave placed on the trajectory computations. 

If, at some future date, it is desired to upgrade the analytical model, 
the inclusion of the * term and the modification of the trajectory 
calculations to three instead of two dimensions is a recommended procedure. 
The total charge exchange ion plume calculation, then, would consist of 
a determination of the neutral density and direction distributions at a 
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given (r,z), including all possible neutral source points on the thruster 
face, and the subsequent trajectory calculations. The final charge 
exchange ion deposition pattern would then require the integral over all 
possible charge exchange locations (r,z). Uncertainties would remain 
in the calculations from the charge exchange ion presence and its 
possible perturbation of the potential. In addition, uncertainties in 
neutral density will be present. These final model elements are discussed 
in Section 5.4.3, which follows. 

5.4.3 Neutral Plume Model Limitations 

The charge exchange ion production calculations utilized a 

neutral plume derived with slmplifyiug assumptions that all neutrals 

emerge from the thruster, that the neutral emission density is uniform 

over the thruster face, and that all atoms emerge with a single velocity, 

^o,th* by the thruster jwall temperature. In addition, the 

angular distribution function for each neutral emission source point has 

2 3 

been chosen as either cos 0, cos 0, or cos 0. An upgrading of the 

analytical model should include possible non-uniformity of neutral 

emission over the thruster face, the inclusion of neutrals released at 

the plasma discharge neutralizer, and the possibility of angular 

2 3 

distributions in emission which depart from cos 0, cos'^0, or cos'^6. A 
complete upgrading of the analytical model is, clearly, a non-trivial 
procedure. 
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6.0 HIGH ENERGY ION FLUX MEASUREMENTS 

6.1 General Considerations 

Measurements of high energy ion fluxes have been obtained from 
the Engine J^, the 1-1/2" the Swinging J_^, and the 4" The high 

energy ion groups of interest here have been previously described in 
Section 1.0 and are both Groups I and II. In the definition used in 
this program, Group I ions are formed in the electron bombardment 
discharge and have been accelerated by the complete potential difference 
between the bombardment plasma and the thrust beam plasma. Group II ions 
are the result of charge transfers between Hg"*" and Hg° in the accelerator 
grid- to-screen grid interspace and at potentials positive with respect 
to the thrust beam plasma. 

The Engine and 1-1/2" probes provide a measurement of 
ion flux within the parabolic core/exponential wing regions. Because of 
cup orientation relative to the axes of rotation by which cup movement 
is attained, neither the Engine or the 1-1/2" probes are effective 
in measuring high divergence angle ion fluxes. For the large divergence 
angles the Swinging and 4" probes are used. The Swinging is 
capable of measurement from the beam axis to 90® divergence angle. 

The 4" has principal regions of effectiveness in the angular divergence 
range above 30®. 

Determinations of energetic ion flux at high angles is not of 
major importance in determining thrust efficiency, since the quantity 
of these particles is small c->.Tipared to the bulk of the (narrowly 
diverging) thrust ions. The major concern raised by the energetic, high 
angle ions is for their possible interception on other system element j 
surfaces or on spacecraft surfaces. In Reference 5 the notion of c = j 
contours is advanced, and, from allowable surface impact and possible ’ 
ion thruster throughput on a primary thrusting mission, location of 

spacecraft surfaces will generally have to be outside of the * 10 cm 

26 ' 
contours. (Note: For 10 ions released by an ion thruster, deposition 

Q O IQ 2 

on the 10”° cm“ contour is 10 ions/ cm and can cause an erosion of 

1000 Angstroms of surface material). In these measurements it will be 

-8 —2 

seen that satisfaction of surface placement at r • 10 cm may not 
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be co-venient for all situation.^, and some procedures to 
angle high energy ion flux will probably be desirable, 
decel will be demonstrated to be an effective method of 
reductions in the large divergence angle ion flux. 

6.2 Engine and 1-1/2" Measurements 


diminish high 
The use of minimum 
introducing 


Measurements of the thrust ion current density at z - 4.7, 

10.0, 15.0, and 20.0 cm are given in Figures 46. 47, 48, and 49. The 
data given there illustrate several commonly observed features in thrust 
ion density. Near the thruster face (z = 4.7 and 10.0 cm) specific ion 
optical features of the discharge chamber and accelerator grid system are 
still evident. Because of these fine structure details, density gradients 
derived from Eq. (10) and Eq. (11) will not be a precise fit to the 
actual density gradients and E fields from the density gradients and 
Eq. (21) will differ from the actual plasma beam internal electric 


fields. 

While the modeled density distribution is not a precise fit, 
it may also be shown to be an adequate representation of the thrust beam. 

The dots in Figures 46 through 49 illustrate calculated values of Eq. (10) 
and (11), the parabolic core/exponential wing density model, with 
Tob = TO cm, a^^ = 5 cm, and = .2. 

Another measured quantity given in Figures 46 through 49 is the 
floating potential in the thrust beam plasma, as determined by electrically 
isolating the Engine elements (across 10 megohms) and recording the 
floating voltage of the probe as it is moved through the beam. There 
are several features of interest in this data. First, the floating 
potential is essentially the known potential on the neutralizer keeper electrode 
indicating good coupling of the neutralizer to the thrust beam. This 
tight coupling and small injection potential indicate, in turn, a 
comparatively low value of T^, thrust beam neutralizing electron 
temperature. This is borne out in the data of Figures 46 - 49 by the 
very weak increments in floating potential as the probe moves from the 
dense plasma near the axis into the more dilute regions of plasma in the 

wing regions. 
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Figure 46. Engine Current and Floating Potential as a Function of 
r at 2 “ 4.7 cm, with Computed Values of Density from Parabolic 
Core/Exponential Wing Thrust Beam Model. 
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Figure 49. Engine J. Current and Floating Potential as a Function 
of r at z ■ 20 cm, with Computed Values of Density from Parabolic 
Core/ Exponential Wing Thrust Beam Model. 



A third aspect of interest in the Engine floating potential 

is the invariance in floating potential for increases in z. Section 5.4.1 
has noted possible failure of to be described by axial density 
gradients in Eq. (10) and the potential formulation of Eq. (21). It 
would appear from the measurements that axial E fields are very weak in 
this region of the beam, and that the charge exchange ion trajectories 
calculated from Eq. (10), (11) and (21), have overestimated the axial 
acceleration of the charge exchange ions. 


The measurements of the 1-1/2" probe are similar to those 
of Engine J^, but allow determinations of larger z values. For 

brevity, results from this probe have not been included here, but are 
given in the Engine Operation Data Volume. 

6.3 High Angle High Energy Ion Measurements 

6.3.1 4" Measurements 

Figures 14 through 23, in Sections 4.5 through 4.7, have 
illustrated the 4" probe signal as a function of z for all ions, for 
ions with energies greater than 95 eV, and for ions with energies less 
than 25 eV, as. engine operation conditions (see Table 3) were varied. The 
discussion in Section 4 concentrated attention on the Group IV ions. 
Discussion in this section will consider the Group I and Group II ions. 

For brevity in this report, the figures will not be repeated in this 
section and reference is made to their earlier presentation. 


Figures 14 and 15 illustrate hard ion currents as the engine 
screen grid potential is varied from 1.5 kV to 0.7 kV. The increase in 
the decel-accel for Engine Condition 18 (Table 3) clearly causes an 
increase in the high angle energetic ion flux. At z » 5 cm (6 -v 80 ) , 
for example, the cup current increases from 4.5 gamperes to 24 gamperes, 
as screen voltage lowers from 1.5 kV to .7 kV, increasing the decel-accel 
ratio from 1.33 to 1.71. (Note that decel-accel ratio. T, is defined here as 
(V^ + jVg|)/Vg and approaches unity as V^ >• 0). 

A second quantity of interest in the data of Figure 14 is the 
magnitude of high angle high energy currents near 0 90“ . Since cup 

area is -- 100 cm^, the 1 gampere of hard ions observed at z 0, correspond 
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to a current density of 10 amperes per square centimeters, and to 

J-) 10~^cm From previous discussion it has been noted 

that maximum allowable e for hard ions, to avoid surface damage, is 
% 10~^cm"^ and, thus, that safe surface placement for this present engine 
could not be in the forward hemisphere. It will be seen that this 
condition can be improved by reductions in accelerator grid potential to 
"minimum" deceleration levels. 

Figures 16 and 17 illustrate hard ion currents as a function of z 

for a nominal 2 ampere thrust beam, 1.1 kV Oi. screen potential, and -.55 and 

-.3 kV as accelerator grid potential. Clearly, the reduction of the 

decel-accel ratio has caused a reduction in high angle energetic ion 

flux. This is demonstrated again in Figures 18, 19, and 20, for a nominal 

1 ampere ion thrust beam, 1.1 kV screen potential, and accelerator grid 

potentials of -.1 kV, -.5 kV, and -.7 kV. At z » 10 cm (6 70*), the 

-7 2-7 

harH ion current density increases from 1 x 10 amperes /cm , to 7 x 10 

9 7 2 

amperes /cm , to 13 x 10 amperes/ cm as decel-accel ratio increases from 
1.09 to 1.45 to 1.64. 

Figures 21 and 22 illustrate hard ion currents as a function of 
z as bombardment discharge potential varies from 43 volts to 34 volts 
(for fixed total discharge power). There is no apparent variation of 
significance in the high angle energetic ion flux as this engine parameter, 
varied, and, on reasonable grounds, none was expected. 

A final measurement of high angle energetic ions by the 4" 
as thruster conditions varied is in Figure 23, where neutralizer flow was 
altered to cause a variation in keeper potential. No significant shift 
in hard ion current was observed, and none was expected. 

6.3.2 Swinging Measurements 

The Swinging probe allows measurement of ion currents at 
divergence angles up to 90*. Because of the mounting of the probe and 
the method of probe motion, the axis of this Faraday cup intersects the 
thruster face at r = z = 0 for all 0 . 

Examples of Swinging total ion current (both hard and soft 

ions) are given in Figures 50 and 51 for thrust ion beams of 1.5 and 


96 






!(/* SWINGING J 




0,5 amperes. (Additional data curves from this probe are in the Engine 
Operation Data). Outside of a central, and comparatively uniform, region 
from Q = 0®to 9 30^ ion flux falls off as v exp(-K0). Two aspects of 

this data should be noted. The first of these is that the probe mounting 
arm length has been set at a comparatively short value so that the probe 
may be rotated to 9 ® 90® without collision with testing chamber walls. 
Because of this short arm length, the finite width of the plasma thrust 
beam appears as an additional angular spread. The true thrust ion 
divergence is, thus, less than that indicated by the 9 = 0° to 0 = 30® 
figure above. A second important aspect to the probe data is that total 
ion current signals are largely dominated by low energy charge exchange 
ions at large divergence angles. The value of K, thus, in the exp {-K0 } 
formulation given above is not readily apparent from total ion current 
measurements, and retarding potential analyses of the probe currents are 
required to determine the hard ion current component. These retarding 
potential analyses will suffer signal-to-noise problems because the hard 
ion currents are minute and occur in the presence of large quantities of 
lower energy ions. 

Figures 52 and 53 illustrate the total ion current and low energy 

and high energy components as a function of 0 for a nominal 2.0 ampere ion 

thrust beam. The total current and soft ion current are given in Figure 52, 

while Figure 53 illustrates the hard ion component. The illustrated case 

is Engine Operatior Condition 19, (see Table 3), a minimum decel condition 

(accelerator grid potential of -.3 kV), which has been shown (by the 

4*’ probe data) to reduce the high angle high energy ion flux. Two 

final examples of Swinging probe data are given in Figures 54 and 55. 

Shown there are the hard ion currents as a function of 9 for 1.0 ampere 

ion thrust beams. Figure 54 is a minimum decel condition (Engine Point 22) 

while Figure 55 illustrates a nominal (V = -.5 kV) decel condition. A 

8 

comparison of the curves demonstrates that the minimum decel condition 
has resulted in a reduction of hard ions at high angles. 

6 , 4 Testing Chamber Ambient Pressure Effects 

Section 4.2 and Figures 9, 10, and 11 have discussed and 
illustrated ambient pressure effects as they tend to alter the measurements 
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Figure 52. Swinging Total Ion Current and Soft Ion Component as a 
Function of 0 for Engine Operation Data Point 19 (Minimum Decel Condition) 
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Figure 55. Swinging J Hard Ion Current as a Function of 0 for 
Engine Operation Data Point 28 (Nominal Decel Condition). 



of low energy ions 4 The presence of the testing chamber gas also causes 
alterations in high energy ion flux at high divergence angles. Figure 56 
illustrates these pressure effects for high angle hard ions. While these 
effects are generally not large, it is possible that some of the remaining 
signals near 6 90® (see, for example Figures 53 and 54) are the results 

of ambient gas which diffuses with the screen grid-to-accelerator grid 
interspace and causes a "pressure effect" Group II ion signal. These 
effects will be discussed further in Section 7. 
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Figure 56* Energetic Ion Current in the 4” Probe as a Function 
of Axial Distance, z, for Two Testing Chamber Conditions for a 
1,0 Ampere Thrust Beam. 




I 


7.0 HIGH ENERGY HIGH ANGLE ION FLUX MODELING 

7 . 1 General Considerations 

Analyses and modeling of high energy ion flux will not be 
concerned with small divergence angle ions. The calculation of ion 
trajectories for the bulk of the Group I ions has received extensive 
treatment elsewhere and will not be discussed further here. The specific 
concern of this section is the hard ion flux for divergence angles 

above 45 , since it is these high divergence particles that create the 

greater part of thruster/spacecraft integration problems. 

Because of the low levels of hard ion signal currents at high 
angles and because of the comparatively high fluxes of Group IV ions in 
these regimes, the retarding potential analyses cannot be carried out over 
the entire range of possible ion energies. Instead, the "hard" ion 
current is defined as those ions possessing in excess of 100 eV and the 
calculations of hard ion flux deposition patterns will be based upon 
ion flux measurements made at % lOO volts retardation. Even though the 
retardation voltage in the measurements has been set at a level which is 
small compared to screen grid potential, it is felt that the ions still 
present in the cup for this retardation setting possess energies, in 
general, significantly above 100 eV. 

7.2 Hard Ion High Angle Flux as a Function of Accel-Decel 

Section 6.3 has discussed the variation of the hard ion flux 
at high angles as screen grid potential and accelerator grid potential 
varied. (See figures in Section 6.3 and also 4.5). The effects are 
clearly evident, particularly for minimum decel conditions. 

It will be advanced here that the bulk of the hard ions seen at 
high angles are Group II ions, have resulted from a charge transfer in 
the screen grid-to-accelerator grid interspace, and have probably been 
subjected to comparatively high decel-accel (that is, their energy as 
they pass through the accelerator grid plane is significantly above the 
final energy the ions possess as they move into the neutralized thrust 
beam region). 
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To examine the comparative regions of accel-decel ratio in the 
interspace between screen and accelerator grids, as a function of 
accelerator and screen grid potentials, the first assumptions are that 

potential in this region is described by space charge limited flow 

■ >“ 

(E = 0 at X * 0, the screen grid location) and that the planar Child- 
Langmuir relation holds. The use of x as the distance variable here, 
instead of the z notation used for axial distance in the parabolic core 
exponential wing model (Eqs. (10) and (11)), is for convenience in the 
sign of the variable and to avoid confusion over the previously chosen 
zero (z = 0 at the accelerator grid). For x = 0 at the screen grid, 

X =s x^ at the accelerator grid, and planar space charge limited ion 
flow in the intervening space, potential is given by 

V(X) = Vg - (V^ + Vg)(^) (26) 


whare is screen grid potential and V is the magnitude of the accelerator 
grid potential. Using the definition of decel-accel ratio as a function 
of X as 


leads to 


V(x) + V 
(x) = S. 


V(x) 


1 - (^ -) 
‘a<") = V ^ 


4/3 


-s (2L_) 

V + V ''X '' 
s g o 


4/3 


(27) 


(28) 


For X = 0, Eq. (28) reduces to the conventional form of decel-accel ratio 
(i.e. (V + V )/V ). The value of x for which T (x) ■> is denoted as 

3/4 


V * X ( ■ ' 

crit o W + V 


(29) 


For O'-^'^'crit’ transfer produces a Group II ion, capable of 

escape into the thrust beam. For . escape is energetically 

forbidden, the ion Is designated as "Group III", and is collected at 
the accelerator grid. 


107 


T 


Using Eq. (26), - 1.1 kV and = .1, .5, and .7 kV, the 

decpl-accel ratio has been computed as a function of x and is illustrated 
in Figure 57. It is apparent from these curves that the decel-accel 
ratio of Group II ions formed in Case 1 (minimum deed condition) is 
significantly less than for Cases 2 and 3. For example, at x = .5x^, 

^ “ 1.160, 2.075, and 2.815 for Cases 1, 2, and 3. Group II ions formed 
at X = * expected, thus, to emerge into the thrust beam 

with successively larger cones of divergence for each of the cases above. 

This does not mean that all Group II ions created at a given x in a given 
(Vg,Vg) case emerge with a single divergence angle. 

For an actual thruster, the Group II formation at x in the 

interval dx, takes place in a cylindrical volume whose radius is 

determined by the radius of the thrust ion beamlet at that point of 

its passage between the screen and accelerator grid. Denoting the outer 

radius of this beamlet as co(x), and noting that the Group II ions formed 

2 

in the cylindrical volume element ttoj dx will be formed at a rate 

^^+cxll 5 

~ “ ( 30 ) 

where J^^^(x) is thrust ion current density at x (and is probably non-uniform 

within the radius w(x)), and >^j^g(x) is Hg° neutral density at x from 

neutrals escaping from the thruster. Using nominal values of x . F . 

2 0 * 0 * 

^o,th* °cx* t ^ total charge exchange formation rate 

of Group II ions of from 10 ^ to a few times lO"^ of I . The total 

+»t 

formation of Group II ions for a 1 ampere thrust beam would range thus 
from a few hundred microamperes to approximately 1 milliampere. Of that 
total production, those ions formed on the axis of the beamlet will 
probably emerge without significant divergence. Group II ions formed at 
the edge of the beamlet, and for which there is a significant departure 
of from unity will be deflected through large angles in the 

resulting passage through the accelerator grid and into the thrust beam. 

It is clearly apparent that regions of significant decel-accel are more 
prevalent in Cases 2 and 3 than in Case 1 and that some of the Group II 
formed at these x values will emerge at high angles. 
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The treatment above is acknowledged to be only an approximate 
one. To improve this description would require the use of the computer 
generated potentials in the regions from the screen grid on into the 
thrust beam, taking Into account thrust ion space charge density and a 
complete description of screen and acceleration grid geometry. In 
addition, more accurately modeled n^^ is required, including, perhaps, 
some inclusion of backstreaming neutrals from the testing chamber. 

Even with these modeled elements complete, the calculation of Group II 
trajectories still appears as difficult and possibly inaccurate because 
of the comparatively reduced ion energy and the consequent severe 
perturbations on the trajectory by all of the electric field patterns 
present in the screen to accelerator grid interspace and in the sheath 
region from the accelerator grid to the thrust beam plasma. A reliance 
upon experimental measurements rather than analyses would appear as 
the most promising approach. 

A final aspect of the Group II calculations to be discussed 

here is the final energy of the Group II ions which have encountered 

significant decel-accel effects and which, presumably, are the Group II 

ions seen at high divergence angles. If it is considered that f >2 

a 

leads to significant divergence, then the maximum Group II ion energy (for high 
divergence angle) in Case 1 is only 100 eV, while in Cases 2 and 3 it 
would be 500 eV and 700 eV. Since the ion sputtering ratio is a 
rapidly rising function of ion energy in this energy regime, the 
sputtering damage per Group II ion at high angles will also rise rapidly 
in moving from the conditions of Case 1 to the conditions of Case 3. 

There are, thus, multiple benef its-~which are obtiiined through the use 
of a minimum decal. Going to minimum decel diminishes the interval 
Ax in x^ for which the higher V values are obtained, and, even within 
those now diminished intervals, has lowered values of the average of 
^ a " * Finally, those Group II ions which do encounter a I' 

a 

sufficiently high to cause severe divergence, emerge into the thrust 
beam with energies of only |V^| . 

^ C alculate d Deposition Contours for Har d Tons 

To assess the possible damage to spacecraft surfaces under hard 
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ion deposition, it has been found convenient to employ c contours, where 

-2 

the units of c are in cm and denotes the flux of a given particle specie 
2 

(in amperes/ cm ) divided by the thrust beam current. Using the notation, 

hard ion deposition per square centimeter during a thrust 

mission in which a total of N. thrust ions are released is e.N. 

+,t + -f,t 

Earlier discussion of allowable deposition have noted that saf^' spacecraft 

—8 “2 

surface placement will generally require £ < 10 cm 

To evaluate the contours for hard ions, the Swinging 
probe data on hard ion flux is used to determine the angular distribution 
function. For convenience in the calculations, it will be assumed that 
all hard ions emerge from r * z = 0. This "point source" approximation 
will be generally valid for most of the relevant regions in r and z for 
spacecraft surface placement. 

Since the Swinging arm length is only several times the engine 
radius, distributed source effects will be present. The use of the Swinging 
flux data and the assumption of point source emission from the thruster 
will result in conservative, upper bound, estimates of the c contour 
placement. 


Examples of contours for hard ions are given in Figures 58, 

59, 60, 61, and 62. Figures 58 and 59 are Engine Operation Conditions 2 

and 19 and are 2 ampere thrust beams under nominal and minimum decel 

conditions, \ftile the use of minimum decel makes surface placement 

somewhat easier by causing a given c contour to move inward toward the 

thrust axis, it should be noted that surface placement outside of the 
*“8 “2 

10 cm contour clearly causes surface placement either at high divergence 

angles near the thruster or at comparatively large axial and radial separation 

distances. Figures 60, 61, and 62 illustrate the contours for a beam of 

1 ampere under accelerator grid voltage conditions of .1, .5, and .7 kV. 

Again, the inward motion of the contours for minimum decel conditions 

is observed and is desirable, but, again, surface placement outside the 
-8 

10 contours requires large angular or spatial separation. 
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Igure 58. Computed Hard Ion Deposition Contours for r, z, and 6 
rom Swinging J Data on Engine Operation Data Point 2 (Nominal Decel). 
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e 60. Computed Hard Ion Deposition Contours for r, z, and 9 
Swinging Data on Engine Operation Data Point 22 (Minimum Decel) 
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gure 61* Computed Hard Ion Deposition Contours for r, z, and 6 
om Swinging J Data on Engine Operation Data Point 23 (Nominal Decel) 
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Figure 62. Computed Hard Ion Deposition Contours for r, z, and 0 from 
Swinging J Data on Engine Operation Data Point 24 (Increased Decel) . 





8.0 EVALUATION OF THRUSTER IN FLIGHT DIAGNOSIS FROM HIGH ANGLE 

ION MEASUREMENTS 

8 . 1 -eneral Considerations 

This section will discuss the use of an ion flux measuring 
probe as an in-flight diagnosis of thruster operation. Several assumptions 
will be made. The first of these is that the total information on thruster 
performance includes measurements of all relevant ion engine voltages and 
currents, plus the outputs of any temperature measuring sensors. The 
ion flux determinations, then, are a complement to a second, and extensive, 
series of ’^terminal** measurements. A second assumption here will be that 
the probe is stationary in position. This assumption is introduced to 
simplify the discussion and to focus attention on diagnostics that are 
of minimum cost. It has been demonstrated in the success of the SERT II 
engine test flight that in-flight measurements of the ion beam flux can 
be carried out with movable Faraday cups and emissive, potential measuring, 
probes. The principal question remaining, then, is whether meaningful 
diagnosis can be accomplished with probes of reduced complexity, which 
leads to the assumed condition of a fixed probe position. 

The assumption of a fixed probe position leads immediately to 

considerations of the erosion of probe surfaces under thrust ion 

interception. From the discussion of previous sections, it has been 

advocated that the placement of surfaces should be outside of the = 

-8 —2 

10 cm contour for energetic ions which leads to a requirement for 
either high angular placement or large radial or small axial separation of 

the probe surfaces from the r = z ® 0 point (at the thruster face center). 

—8 —2 

This advocacy of location outside the 10 cm for energetic ions may 
be overly conservative, since the high angle ions are probably Group 11 
rather than Group 1, and, because of reduced energy in Group 11 compared 
to the thrust ions, will have lowered values of ion sputtering. 

A final aspect to be treated in these general considerations 
is the focus of the probe diagnosis. Although the probe to be advocated 
for spacecraft use will possess the capability for several different 
measurements, the principal Interest in the diagnosis will be for Group IV 
ions, from which engine performance in propellant utilization will be 
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deduced. The secondary target for the probe diagnosis will be to determine 
the effectiveness of beam neutralization, and the tertiary goal will be 
to determine the high angle hard ion flux. This final measurement may 
be of use in evaluating the in-flight wear-in of the thruster grids. 
Conclusions as to wear-in effects and hard ion count are only tentative, 
since the present measurements program has not been directed toward long 
term alteratxons in the shapes and magnitudes of the several species of 
ion plumes from the thruster, 

8. 2 Probe Placement 

Probe placement is principally determined by the location of the 
peak in the genuine Group IV ion flux. From Figure 12 and Figures 14 
through 23 (and for a probe similar to the 4" see Section 8,3, which 
follows), the optimum location in z (for r 28 cm) is in the range 
from 10 cm to 20 cm. For probe location outside this range, 
increasing problems will be encountered between laboratory and in-flight 
results since the regimes of z < 10 cm and z > 20 cm are increasingly 
dominated by facility effect Group IV ions. Another reason for probe 
placement in this ^’plateau** region is that this region, because of its 
"flatness" (small variations in probe current for variations in z 
location) will probably not be subject to probe misinterpretation because 
of minor changes in the shape of the Group IV ion plume. For locations 
of the probe near z = 10 cm or for z < 10 cm, the rapid drop-off in 
probe signal for motion in z toward z = 0, leads to a condition of 
increased sensitivity to plume shape changes as contrasted to the desired 
measurements of plume magnitude , 

A major concern in probe placement is the extent of the energetic 
ion flux. The erosion of probe surfaces (and possible secondary mass 
transport and deposition effects) clearly weighs against probe placement 
at z - 20 cm. It is also not desirable for the Group IV determination 
to be carried out in the presence of comparatively large Group I and II 
ion fluxes. An acceptable condition can be equal levels of energetic 
and soft ion fluxes (these can be accurately separated by retarding 
potential analyses). These equal levels of hard and soft ions arc 
encountered near z = 10 cm for nominal accel-decel conditions, and near 
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z = 20 cm for minimum decel conditions. The use of a given accel-decel 
ratio may depend upon many different thruster and mission considerations, 
and it does not appear advisable to allow plume diagnosis requirements 
to be the "driver" in this area. Plume diagnosis requirements should be 
considered in the selection of an accel-decel condition, and, in a 
self-consistent manner, probe placement will depend upon the selected 
accel-decel condition. Depending upon many factors which may see some 
future alteration, it would presently appear that probe placement should 
be in the range 10 cm < z <20 cm (for r 28 cm). Section 8.3, which 
follows, will consider probe orientation at the selected location, and 
probe configurations and capabilities. 

8*3 Probe Configuration 

Section 8.1 has taken the position that the principal focus of 
the probe measurements will be upon the magnitude of the charge exchange 
ion flux, from which (using Eq. (18) and appropriate calibration 
constants) can determine propellant utilization (essentially term F 
in Eq. (18)). For these Group IV flux determinations to be as insensitive 
as possible to other effects, such as small alterations in plume shape 
or in charge exchange Ion trajectory, it is essential that the probe 
be a "total flux" measuring probe, that is, that the probe construction 
and orientation be such as to have a large solid angle for charge exchange 
ion trajectory acceptance. For this reason it is recommended that the 
probe have a conf igurction similar to the 4" probe, with the orientation 
of the surface normal of the collector such that it passes through the 
r = z = 0 point at the center of the thruster face. 

Since both hard and soft ions will be directed into the cup and 
since determinations of both ion species are of interest, the probe 
should be a raulti-grldded probe (at least two grids are required) capable 
of carrying out a retarding potential analysis of the Incoming ion flux. 
This requires an electronics package capable of applying a varying 
(positive) bias on the retarding potential grid from 0 volts to, at 
least, +100 volts. The outer grid must be negatively biased to prevent 
drainage of electrons from the plasma plume to the inner, retarding 
potential, grid. A third grid may also be useful as a means of 
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suppressing secondary electrons and photoelectrons from the collector 
surface. This third grid, between the retarding potential grid and the 
collector surface, would be biased negatively with respect to spacecraft 
ground. The collector surface would be, essentially, at spacecraft ground 
potential. 

A final useful output of the probe (in addition to soft and hard 
ion fluxes) is plasma thrust beam potential. For a probe of the size 
recommended, and for the expected flux magnitudes, the probe can be used 
effectively for floating potential measurements with only moderate values 
of isolation impedance. The floating potential measurements would be used, 
in turn, to determine in-flight effectiveness of the thrust beam neutralizer. 
In the floating potential mode of operation, all probe components (case, 
outer grid, RPA grid, inner grid, and collector) are connected together 
and are electrically isolated from spacecraft ground by an isolation 
resistance of at least 10 megohms. 

Table 5 presents recommended probe configuration, size, and 

placement. 

Table 5. Configuration Details, Size, and 
Placement of Ion Thruster Plume Diagnosis Probe. 


Probe Location 

Probe Orientation 

Collector Area 
Grids 

Outer Grid Potential 

Retarding Grid Potential 

Suppressor Grid Potential 

Collector Measurement Impedance 

Configuration for Floating 
Potential Determinations 


10 cm < z 20 cm 

r 30 cm 

Collector surface normal 
passes through r = z = 0 

50 - 100 cm^ 

3 

-10 to -20 volts 
(fixed) 

. 0 to +100 volts 
(varial)le in steps) 

-10 volts 
(fixed) 

10^ : to 10^ 

(variable in decade steps) 

All probe elements (conm;ct:e(l 
together and isolated f rom 
spacecraft ground hv at least 10 
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8.4 


Multiple Thruster (Cluster) Effects 


The measurements program carried out in the 5' x 11' facility 
has determined the plume characteristics of only a single ion engine. 

From Section 5.4, it should be expected that the presence of a thruster 
and its charge exchange plasma plume can be affected by the presence of 
a second thruster, and accompanying plasma plume. The extent and 
particular characteristics of these pile-up effects on the Group IV ion 
trajectories are not known at present and should be an element in future 
beam diagnostics programs. Because of facility effect ion production, 
and because the operation of multiple thrusters will create additional 
pumping problems for the test facility, these "cluster" effect measurements 
should be carried out J.n the largest and highest pumping speed facilities 
possible. 
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SUMMARY AND RECOMMENDATIONS 


The ion flux from a 30 cm ion thruster has been examined for 
Group I, Group II, and Group IV ion species. The presence of the testing 
facility boundaries has no significant effect on Group 1 and Group II ions. 

The Group IV ion measurements, however, are subject to influence by 
facility effects. By reducing chamber pressure into the range below 
4 uTorr, the principal features of the charge exchange ion plume may be 
determined. Facility effects are still present at these chamber pressures, 
however, and even further reductions in chamber pressure are recommended. 
Depending upon point of observation, direction of observation, or range 
of ion energy, measurements may be influenced by the ambient Hg® density, 
even for the highest pumping speeds and the lowest possible chamber 
pressures. Variation of facility pressure and correlation of facility 
pressure with observed plume fluxes can be used as an aid to identifying 
possible facility effects. This method is not precise, however, unless 
an accurate ^ situ determination of all partial pressure contributions 
is made at the given test facility total pressure. 

The Hg"^ charge exchange ion plume has been investigated under 

a series of engine operation conditions. For thrusters with discharge 

conditions set for efficient ionization, the Group IV flux scales as the 

thrust ion current, I . These results necessarily imply an improvement 

+ , t 

in propellant utilization as thrust ion current is increased. Since 
propellant utilization efficiency varies as thrust ion current varies, 
these experiments determine the results of an explicit variation in ^ 
and an implicit variation in propellant utilization. The charge exchange 
production rate was also examined as an explicit function of propellant 
utilization by holding I fixed and varying the emission current in the 
bombardment discharge (which varies propellant utilization) . Charge 
exchange Ion flux increased for diminished propellant utilization, as 
expected from the charge exchange ion production model. The experiments, 
however, did not examine deliberately ’’spoiled" utilization conditions 
(defined here as a significant lowering of propellant utilization away from 
the values obtained for normal engine operation at a given thrust inn 
current). Since a major possible use of Group IV measuremimt s in space 
should be to in-flight monitoring of propellant utilization, it is recommended 
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that additional Group IV moasuromcnts be made under bombardment discharBc 
conditions leading to deliberately "spoiled" utili..ation efficiency. 

The charge exchange ion plume shape has not exhibited any major 
variations as screen grid potentials and accelerator grid potentials move 
through relevant ranges. Changes in discharge bombardment voltage (for 
constant discharge power) have not revealed variations in Group IV plume 
shape. Both shape and magnitude of the Group IV plume were invariant 
to alterations in neutralizer heater condition (within the range of 
variation utilized) . 

The results above appear to confirm analytical models of charge 
exchange ion production and deposition. In the computed charge exchange 
ion trajectories, the governing situation is the shape of the major 
portion of Group I ions. If this modeled beam expands, the Group IV 
plume appears to have a counter motion, while a narrowing Group 1 beam 
leads to an expanded Group IV plume. 

The total deposition pattern of these Group IV ions has not been 
computed. Such computations are possible, but will require additional 
specification of model parameters if results are to be precis 


The Group I and Group II ions 
divergence angle range from 0 to 90 . 


Viave been examined in the 

The cutoff in Group II abundance 


near 90° is more abrupt for reduced deceleration, and several Important 


reductions in hard ion deposition effects can be obtained through the 
use of a minimum decel condition. These factors should be included with 
other engine operation data in the ultimate choice of thruster operation 

parameters. 


The recommended properties for an in-flight diagnostic probe 
have been given (Table 5) and include retarding potential analysis 
capability, electrical floating capability (for measurements of potential 
in the neutralized thrust beam), and a broad range of acceptance angles 
for ion trajectories. The location of this probe should be at high 
angles (■ 60° to • 80°), with emphasis on Group TV Hg"^ measuremntits , 
but with a capability, also, for Group 11 identification. 
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The measurements of beam efflux have taken place over only a 
limited period of thruster operation so that there is, as yet, no 
firm evidence for the plume shapes and magnitudes for ultimate, run-in 
conditions. It is recommended that some of the long duration thruster 
test runs have these diagnostic tests performed during the total period 
of beam release. It is also recommended that plume measurements for 
single thrusters and groups of thrusters be carried out in the largest 
size and highest pumping speed chambers available to more clearly define 
the cluster effects and to permit determinations, over broader regions 
of parameter space, of the genuine engine effluxes. 
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